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Mr.  C.  E.  Chatham,  Chief,  NHD;  and  under  the  general  supervision  of  Mr.  Charles 
C.  Calhoun,  Jr.,  Assistant  Director,  CHL,  and  Dr.  James  R.  Houston,  Director, 

CHL.  Word  Processing  and  formatting  were  done  by  Ms.  Myra  E.  Willis,  CHL. 
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Summary 


Introduction 


The  U.S.  Army  Engineer  Waterways  Experiment  Station  (WES),  Coastal  and 
Hydraulics  Laboratory,  provided  technical  support  to  the  U.S.  Army  Engineer  Divi¬ 
sion,  Pacific  Ocean  (POD)  for  the  study  of  Kaumalapau  Harbor.  Kaumalapau  is  a 
privately  owned  small  barge  harbor  located  on  the  southwest  coast  of  Lanai,  Hawaii, 
and  is  the  only  commercial  harbor  serving  the  island.  Fuel  and  commodities  for  the 
island’s  inhabitants  come  through  this  harbor.  The  harbor  is  semicircular  in  shape 
and  has  no  distinct  entrance  channel.  A  wharf  is  located  along  the  northern  portion 
of  the  harbor  and  is  protected  by  a  76-m-  (250-ft-)  long  breakwater.  The  rubble- 
mound  breakwater  has  been  damaged  by  past  hurricanes  and  provides  minimal 
shelter  to  the  wharf  area. 

Four  interrelated  aspects  of  the  study  were  conducted  by  WES;  (a)  review  of 
wave  hindcast  data,  (b)  field  wave  gauging,  (c)  numerical  model  simulations,  and 
(d)  physical  model  studies. 

Project  Goals  and  Restrictions 

The  goal  of  the  study  was  to  decrease  wave  heights  along  the  wharf  area  without 
adversely  impacting  harbor  navigation.  Initially,  harbor  pilots  were  contacted  to 
ascertain  their  navigational  needs.  Although  the  usual  criteria  established  for 
harbors  is  for  the  entrance  channel  to  be  three  to  ten  times  the  boat  width,  at 
Kaumalapau,  the  vessel  has  already  started  to  make  its  turn  toward  the  wharf  by  the 
time  it  reaches  the  harbor  breakwater  and  thus  the  standard  convention  for  channel 
width  is  not  applicable.  The  pilots  stated  that  the  harbor  is  so  small  and  there  is 
little  room  inside  the  harbor  to  turn  so  that  no  protective  structure  should  encroach 
on  the  navigation  channel.  This  criteria  eliminated  most  protective  structures, 
except  for  an  extension  to  the  existing  breakwater  of  approximately  46  m  (150  ft). 

POD  wanted  to  improve  wave  conditions  at  Kaumalapau  Harbor  so  that  the 
harbor’s  use  would  not  be  limited  by  wave  conditions.  The  established  project  goal 
with  end-user  input,  was  to  decrease  the  incident  wave  height  of  4.6  m  (15  ft)  to 
1.5  m  (5  ft)  at  the  wharf. 


Wave  Climate 


Wave  conditions  and  southern  swell  information  needed  for  numerical  and  physi¬ 
cal  model  input  were  queried  from  Wave  Information  Study  (WIS)  results.  WIS 
analyzed  20  years  of  meteorological  data  to  predict  (hindcast)  wave  heights  every 
3  hr  for  20  years  by  running  a  computer  program  that  infers  wind  speed  and  its 
impact  on  wave  generation  in  the  open  ocean.  In  addition  to  hindcast  wave  data, 
measured  wind  and  wave  data  at  three  nondirectional  National  Oceanic  and  Atmo¬ 
spheric  Administration  (NOAA)  buoys  south  and  southwest  of  Kaumalapau  were 
analyzed. 

WIS  data  and  NOAA  buoys  were  located  a  sufficient  distance  from  Kaumalapau 
Harbor  so  that  wave  data  at  the  site  were  needed  to  calibrate  the  numerical  model. 
Validation  of  a  harbor  model  involves  driving  the  model  with  measured  input  data, 
and  comparing  the  model's  output  to  measured  response  at  one  or  more  sites  within 
the  harbor.  Two  wave  gauges  were  required  based  on  preliminary  results  from  the 
numerical  model  HARBD  and  logistic  constraints.  Data  from  the  field  measure¬ 
ments  were  used  to  adjust  variables  used  in  the  numerical  model,  thus  providing  a 
calibrated  and  validated  model. 

Numerical  Model  Studies 

The  numerical  model  studies  had  four  main  objectives:  (a)  assist  in  optimizing 
placement  of  field  wave  gauges;  (b)  validate  the  model  with  field  data;  (c)  maximize 
efficiency  of  physical  model  experiments  by  identifying  most  damaging  incident 
wave  directions  and  most  promising  harbor  modification  plans;  and  (d)  evaluate 
harbor  resonance  characteristics  of  the  final  WES-recommended  plan  relative  to  the 
existing  harbor. 

To  perform  the  wave  response  study  for  Kaumalapau  Harbor,  WES's  numerical 
wave  response  model  HARBD  was  used  (Chen  1986,  Chen  and  Houston  1987, 
Lillycrop  et  al.  1993,  Thompson  and  Hadley  1995).  HARBD  is  used  in  the  calcula¬ 
tion  of  linear  wave  response  in  harbors  of  arbitrary  shape  and  depth.  The  model 
includes  capabilities  for  modeling  both  long  and  short  waves,  bottom  friction,  and 
partially  reflective  boundaries.  The  model  does  not  include  wave  transmission, 
wave  overtopping,  currents,  wave  breaking,  or  nonlinear  effects.  Diffraction  around 
structure  ends  is  represented  as  diffraction  around  a  blunt  vertical  wall  with  a  speci¬ 
fied  reflection  coefficient.  Despite  these  limitations,  HARBD  still  was  suitable  for 
meeting  the  numerical  modeling  objectives  of  the  Kaumalapau  Harbor  study. 

HARBD's  principal  output  consists  of  an  amplification  factor  and  phase. 
Amplification  factors  are  easily  interpreted  and  predict  those  wave  periods  that 
cause  oscillation  or  tilting  of  the  free  surface  which  in  turn  make  mooring  of  vessels 
difficult.  Phase  information  is  useful  in  long  wave  studies  for  determining  relative 
phase  differences  within  the  harbor  and  interpreting  harbor  oscillation  patterns  and 
which  mooring  locations  to  avoid. 


HARBD  studies  were  conducted  to  compare  the  existing  and  recommended 
harbor  plan  as  well  as  evaluating  alternative  designs  for  wind  waves  (6-22  sec) 
and  long-period  waves  (23-500  sec).  The  harbor  response  to  long-period  waves 
was  simulated  to  identify  resonant  nodes  of  free  surface  oscillation.  Compared  to 
the  existing  plan,  resonant  nodes  at  several  locations  along  the  wharf  for  the 
recommended  plan  show  fewer,  but  generally  higher,  oscillations  of  the  free 
surface.  The  model  revealed  that  the  shorter  wind  waves  are  smaller  along  the 
wharf  area  for  the  recommended  plan  in  comparison  with  existing  conditions. 

Physical  Model  Studies 

An  undistorted,  three-dimensional  physical  model  of  Kaumalapau  Harbor  was 
constructed  to  further  evaluate  harbor  response  to  short-period  waves  for  existing 
conditions  and  two  alternative  harbor  layout  configurations  for  three  incident  wave 
directions.  Once  the  recommended  configuration  was  defined,  physical  model 
experiments  were  conducted  to  optimize  breakwater  stability. 

Wave  response  at  the  pier  area  was  analyzed  for  three  breakwater  configura¬ 
tions.  Existing  breakwater  conditions  were  evaluated  first  to  provide  a  baseline 
against  which  the  alternative  configurations  could  be  compared.  The  second 
breakwater  configuration  studied  was  the  pre-Hurricane  Ewa  breakwater  length  of 
122  m,  which  extends  to  the  edge  of  the  navigation  channel  (an  additional  50  m 
past  existing  conditions).  The  last  alternative  studied  is  also  the  preferred  plan 
and  was  termed  the  dogleg  breakwater.  The  last  15  m  of  the  breakwater  center 
line  is  aligned  30  deg  toward  the  inside  of  the  harbor  as  measured  from  the 
center  line  of  the  straight  breakwater. 

Of  the  alternatives  studied  and  compared  to  existing  conditions,  the  dogleg 
breakwater  provides  better  attenuation  of  short-period  wave  energy.  The  largest 
waves  outside  the  harbor  that  will  meet  the  1.5-m  (5-ft)  criterion  along  the  wharf 
for  the  existing,  straight,  and  dogleg  breakwaters  (for  waves  from  221  and 
251  deg  measured  clockwise  from  TN)  are  2.0,  2.1,  and  2.3  m  (6.5,  7.0,  and 
7.6  ft),  respectively,  when  averaged  over  the  two  directions  and  numerous  wave 
heights  and  periods.  This  does  not  meet  the  original  design  criteria  established  by 
the  state,  but  these  breakwaters  do  meet  the  criteria  established  by  the  end  users. 
Further  wave  reduction  could  be  met,  but  the  breakwater  would  not  meet  the 
navigational  criteria  and  economic  constraints. 

Increasing  the  breakwater  length  provides  more  protection,  but  care  is  needed 
to  build  the  structure  with  proper  armor  units  to  ensure  durability  and  lower  main¬ 
tenance  costs.  Armor  units  placed  on  the  ends  of  breakwaters  often  are  vulnera¬ 
ble  to  waves  which  overtop  and  diffract  around  the  structure  because  the  units  are 
subjected  to  forces  from  the  opposing  direction  of  a  two-dimensional  situation. 
Heavier  units  often  are  required  for  breakwater  stability  on  the  head  (seaward 
end)  of  the  structure. 

The  design  wave  was  specified  in  a  reconnaissance  report  as  a  9.8-sec,  8.5-m 
(27.9-ft)  deepwater  wave  (U.S.  Army  Engineer  District,  Honolulu  1993)  based 
upon  the  “worst  case”  hurricane  passing  the  island.  To  determine  a  factor  of 
safety  for  the  structure,  wave  heights  exceeding  the  design  height  were  generated 


for  9.8-  and  12-sec  periods.  Wave  heights  associated  with  16-sec  periods  from 
the  221-  and  251-deg  directions  are  smaller  in  magnitude  (limitations  of  the  wave 
machine);  therefore,  only  heights  up  to  Ha  =  8.5  m  were  generated  for  16-sec 
waves.  All  wave  conditions  were  reproduced  at  a  storm  surge  water  level  of 
+  1.52  m  (4.98  ft)  mean  lower  low  water  and  a  duration  of  24  hr. 

Experiment  results  revealed  that:  a  toe  trench  at  the  shoreward  end  on  the  sea 
side  of  the  breakwater  was  necessary  to  obtain  stability  and  densely  placed 
18,144-kg  CORE-LOCS™  (hereafter  referred  to  as  core-loc)  were  stable  for  de¬ 
sign  waves,  while  heavier  core-locs  (31,389-kg)  were  required  for  wave  heights 
exceeding  the  design  condition. 
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Background 


At  the  request  of  the  U.S.  Army  Corps  of  Engineers,  Pacific  Ocean  Division 
(POD),  the  U.S.  Army  Engineer  Waterways  Experiment  Station's  (WES)  Coastal 
and  Hydraulics  Laboratory  carried  out  a  wave  response  study  for  Kaumalapau 
Harbor,  Lanai,  Hawaii.  Kaumalapau  Harbor,  Lanai's  main  commercial  harbor,  is 
being  studied  for  development.  It  is  protected  from  the  open  ocean  by  a  single 
rubble-mound  breakwater. 

The  breakwater,  particularly  the  seaward  end,  has  suffered  damage  from 
previous  storms.  Most  recent  damage  was  caused  by  Hurricane  Iniki  and  winter 
storm  waves  in  February  1993.  In  its  present  configuration,  the  most  seaward 
22.9  m  (75  ft)  of  breakwater  length  has  been  reduced  to  a  submerged  mound.  The 
mound  extends  south  and  slightly  seaward  from  the  existing  breakwater  head.  The 
shallowest  depth  over  the  submerged  breakwater  head  is  approximately  3.6  m 
(12  ft). 

There  is  concern  that  die  existing  breakwater  does  not  provide  adequate  protec¬ 
tion  of  the  harbor  from  approaching  ocean  waves.  Wave  action  in  the  vicinity  of  the 
main  dock  is  the  primary  concern.  Difficult  wave  conditions  are  reported  to  occur 
during  the  winter  season,  particularly  during  high  energy  swell  from  the  north. 


Study  Location 

Kaumalapau  Harbor  is  located  in  a  small  embayment  along  the  south  central  part 
of  the  west  coast  of  the  island  of  Lanai  (Figure  1).  The  general  exposure  is  to  the 
west.  The  northwest  lobe  of  Lanai  shelters  the  site  from  the  north  and  north- 
northwest.  The  entrance  to  Kaumalapau  Harbor  is  formed  by  a  rocky  point  on  the 
soudi  side  and  the  main  Kaumalapau  breakwater  on  the  north  side  (Figure  2). 
Navigation  lights  are  located  on  the  breakwater  tip  and  the  point  south  of  the  harbor. 
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Figure  1 .  Study  location 


The  harbor  entrance  is  a  183-m-  (600-ft-)  wide  opening  into  a  semi-protected, 
40,500-m2  (10-acre)  berthing  area.  Water  depths  in  the  berthing  area  range  from 
9. 1  to  15.2  m  (30  to  50  ft).  The  rubble-mound  breakwater  is  approximately  76.2  m 
(250  ft)  long  with  a  crest  elevation  of  about  3.0  m  (10  ft).  Commercial  operations 
occur  along  a  single  121.9-m  (400-ft)  wharf  in  the  lee  of  the  breakwater.  The  wharf 
is  concrete  atop  a  rock  base.  The  wharf  face  is  concrete  down  to  about  the  waterline 
and  pile  supported  under  water.  Behind  the  piles  is  a  nearly  vertical  rock  face. 

Kaumalapau  Harbor  is  Lanai's  main  commercial  harbor.  Lanai's  second  harbor 
of  note  is  a  small  boat  harbor  at  Manele  Bay,  on  the  southeast  coast.  Primary 
cargos  at  Kaumalapau  have  traditionally  been  outbound  pineapples  and  inbound 
supplies  needed  for  people  and  operations  involved  in  pineapple  farming.  Pineapple 
farming  on  Lanai  has  virtually  ended,  and  facilities  for  recreation  and  tourism  are 
developing.  Kaumalapau  Harbor  will  continue  to  be  the  critical  sea  link  between 
Lanai's  residents  and  visitors  and  the  rest  of  the  world. 

Most  of  the  shoreline  of  Kaumalapau  Harbor  consists  of  rocky  cliffs.  Water 
depth  over  the  harbor  is  generally  between  6. 1  and  15.2  m  (20  and  50  ft).  A  stream 
discharges  into  Kaumalapau  Harbor  in  the  center  of  the  eastern  shore  of  the 
embayment.  The  bottom  is  shallow  and  gently  sloped  in  the  vicinity  of  the  stream 
outlet.  Bottom  material  in  the  harbor  is  generally  sand  with  some  rocky  features 
near  the  shore  and  in  the  east-central  part  of  the  embayment. 
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Figure  2.  Kaumalapau  Harbor 


Harbor  Operational  Considerations 


Before  the  study  commenced,  the  harbor  end  users  -  Sause  Brothers  (fuel  barge 
operators)  and  Young  Brothers  (container  ship  owners)  were  contacted  to  ascertain 
their  operational  and  navigational  needs.  Kaumalapau  is  a  small  harbor  and  the  first 
stipulation  imposed  by  the  harbor  pilots  was  that  no  protective  structure  should 
encroach  into  the  existing  navigation  channel.  Kaumalapau  is  unique  in  that  the 
embayment  is  very  narrow.  The  distance  from  the  6. 1-m  (20-ft)  depth  contour  along 
the  shoreline  adjacent  to  the  cliff  walls  to  the  20-m  (65.6-ft)  depth  contour  of  the 
existing  breakwater  is  approximately  150  m  (500  ft).  The  navigation  channel  is 
indicated  by  the  navigation  buoys  shown  in  Figure  2. 

The  largest  vessel  calling  at  Kaumalapau  Harbor  is  a  fuel  barge  with  a  length  of 
83  m  (272  ft),  a  beam  width  of  23  m  (76  ft),  and  a  draft  of  5.3  m  (17.4  ft).  The  fuel 
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barge  is  accompanied  by  two  tugs,  which  are  30.5  m  (100  ft)  in  length.  The  tugs 
operate  on  the  same  side  of  the  barge  and  the  combined  width  of  the  three  vessels  is 
approximately  30.5  m  (100  ft). 

Container  barges  call  on  the  harbor  on  a  regular  basis.  These  vessels  are  69  m 
(226  ft)  long,  with  a  beam  of  18  m  (58  ft),  and  a  draft  of  4  m  (13  ft).  Sause 
Brothers  has  indicated  that  the  fuel  barge  requires  smaller  wave  conditions  than  the 
container  barge  for  safe  navigation  and  mooring. 

The  operational  condition  established  by  POD,  after  discussion  with  the  end 
users,  was  for  1.5-m  (4.9-ft)  wave  heights  or  less  at  the  wharf  when  4.6-m  (15.1-ft) 
waves  are  present  outside  the  harbor.  These  waves  were  not  linked  to  a  specific 
direction  or  wave  period. 
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2  Wave  Climate1 


A  large-scale  physical  model  of  the  harbor  was  constructed  by  WES  to  opti¬ 
mize  needed  modifications  to  the  existing  harbor.  Wave  conditions  required  to 
drive  the  model  were  requested  from  the  Wave  Information  Study  (WIS). 
Measured  wind  and  wave  data  at  three  nondirectional  National  Oceanic  and 
Atmospheric  Administration  (NOAA)  buoys  south  and  southwest  of  Kaumalapau 
were  analyzed  to  supplement  WIS  data. 

Statistical  products  (percent  occurrence  tables  and  wave  rose  diagrams) 
describing  hindcast  results  at  WIS  stations  34  and  35  were  extracted  from  the  WIS 
Pacific  Phase  I  final  report  (1986),  and  guidance  for  the  use  of  the  products  was 
provided.  A  similar  statistical  analysis  was  performed  with  measured  data  from 
NOAA  buoys  51002,  51003,  and  51005  (Figure  3).  Because  these  buoys  do  not 
measure  wave  direction,  recorded  wind  direction  data  were  used  to  estimate  wave 
direction. 

Results  from  two  Pacific  Phase  I  hindcast  stations  indicated  that,  for  the 
20  years  hindcast,  the  mean  wave  height  is  near  2.5  m  (8.2  ft).  The  largest 
significant  wave  height  is  6.7  m  (22.0  ft),  and  the  mean  peak  period  is  9.2  sec.  A 
total  of  58,440  cases  were  analyzed,  and  99  percent  of  the  waves  came  from  the 
north-northwest  through  north  to  north-northeast,  meaning  that  only  1  percent  or 
less  came  from  directions  that  would  affect  Kaumalapau.  This  result  is  explained 
by  the  shape  of  the  WIS  Pacific  Phase  I  grid,  which  allows  only  a  small  amount  of 
energy  from  the  south  and  west  of  Hawaii  to  come  into  the  hindcast  grid. 

Analysis  of  measurements  from  three  nondirectional  NOAA  buoys  near  Hawaii 
supplements  WIS  hindcast  data.  Occurrence  of  waves  from  south  and  west  of 
Kaumalapau  is  estimated  using  observations  from  a  wind  direction  window  (180- 
to  315-deg  azimuth).  Occurrences  of  total  observations  for  buoys  51002,  51003, 
and  51005  from  that  window  are  1.6,  3.7,  and  7.0  percent,  respectively.  Summer 
swell  (only  0. 1  or  0.2  percent  of  total  observations)  is  isolated  by  limiting  the 
occurrences  to  April  through  September  and  wave  period  to  12.5  sec  or  more. 
Total  wave  observations  are  classified  by  height  and  period,  as  is  swell  for  all 
seasons  and  summer  swell  for  all  directions.  Wave  height  maximums  are  6.2  m 
(20.3  ft)  for  buoy  51002,  6.8  m  (22.3  ft)  for  buoy  51003,  and  5.4  m  (17.7  ft)  for 
buoy  51005  for  the  individual  period  of  record. 
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Field  Wave  Measurements1 


Planning 

Wave  measurements  were  made  at  Kaumalapau  Harbor  to  provide  data  to 
validate  numerical  (HARBD)  and  physical  models  of  the  harbor  response  to  incident 
wave  energy.  (A  secondary  objective  was  validation  of  the  structure  stability  model. 
Since  it  would  require  measurements  of  an  event  of  sufficient  magnitude  to  cause 
quantifiable  damage  to  the  structure,  this  objective  was  not  expected  to  be  met 
during  a  limited-duration  study.)  Validation  of  a  harbor  model  involves  driving  the 
model  with  measured  input  data,  and  comparing  the  model's  output  to  the  measured 
response  at  one  or  more  sites  within  the  harbor.  Planning  the  measurement  program 
requires  specifying  the  location,  duration,  and  type  of  data  collected.  Ideally, 
incident  measurements  coincide  with  the  outer  boundary  of  the  model,  and  there  are 
sufficient  interior  measurements  to  define  the  spatial  variability  within  the  harbor. 

It  is  desirable  to  continue  measurements  long  enough  to  obtain  a  broad  range  of 
incident  conditions  -  up  to  or  exceeding  design  conditions.  Finally,  the  types  of  data 
(wave  energy,  wave  direction,  currents,  etc.)  and  the  range  of  frequencies  measured 
will  equal  or  exceed  the  requirements  of  the  model.  Practically,  fiscal,  logistic,  or 
schedule  limits  will  constrain  the  measurement  plan. 

A  minimum  of  two  sites  was  required;  one  for  incident  conditions,  and  one  for 
harbor  response.  Wave  gauge  placement  was  based  on  preliminary  (uncalibrated) 
results  from  HARBD  and  logistic  constraints.  Both  the  numerical  and  physical 
models  include  water  level,  incident  wave  energy,  frequency,  and  direction  as  input, 
and  the  preliminary  HARBD  results  showed  the  harbor  response  was  sensitive  to 
incident  wave  direction.  The  available  budget  only  permitted  installation  of  non- 
directional  wave  gauges.  Funds  from  other  programs  were  anticipated  to  support  a 
directional  wave  gauge  in  deep  water,  west  of  Lanai,  at  some  point  during  the  study. 

While  depths  at  outer  boundaries  of  the  models  ranged  from  about  29  to  45.7  m 
(95. 1  to  150  ft)  mean  lower  low  water  (mllw),  instrumental  (pressure  attenuation 
factor)  and  installation  (divers'  bottom  times)  constraints  limited  gauge  depths  to  a 
maximum  of  about  2 1  m.  A  site  between  the  toe  of  the  breakwater  head  and  the 
entrance  channel  met  the  depth  constraint  and  provided  some  protection  from  the 
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risk  of  damage  from  tow  bridles.  However,  measurements  at  this  site  will  contain 
some  reflected  energy;  in  fact,  the  HARBD  results  showed  the  entire  model  domain 
was  affected  by  reflected  energy  from  the  breakwater  and/or  the  interior  of  the 
harbor.  Because  of  the  steep  shoreline,  an  "incident  only"  wave  measurement  would 
require  a  gauge  on  the  order  of  miles  offshore.  For  the  single  interior  gauge,  a 
position  on  the  wharf  is  obviously  of  most  interest  from  an  operations  viewpoint. 
The  middle  to  western  section  of  the  wharf  was  recommended  by  the  HARBD 
results,  since  it  showed  unexpected  amplification  for  the  longer  wind  waves  (swell), 
and  contained  fewer  nodes  for  harbor  oscillation  than  the  eastern  end. 

System  Design 


Potential  damage  from  large  waves,  vessel  traffic,  and  vandalism  made  bottom- 
mounted  pressure  sensors  the  best  sensor  choice.  Oscillating  quartz  transducers 
(ParosScientific)  were  used  in  PVC  pressure  housings.  The  transducers  were 
sampled  at  1  Hz,  which  is  sufficient  to  permit  measurement  of  wind-wave  energy  up 
to  0.5  Hz.  While  continuous  sampling  is  desirable,  limitations  of  memory  and 
battery  power  prescribe  a  discrete  sampling  scheme.  Both  sensors  were  sampled  for 
2,048  sec  every  hour.  When  conditions  were  mild  (significant  wave  height  below 
1.5  m  (4.9  ft))  sampling  was  reduced  to  4-hr  intervals. 

Two  data  capture  modes  were  used:  telemetry  to  WES,  to  monitor  system  per¬ 
formance  and  ensure  recovery  of  transmitted  data  in  the  event  of  catastrophic 
damage  to  the  system;  and  internal  storage,  to  retain  data  for  later  recovery  in  the 
event  of  loss  of  the  telemetry  link.  The  telemetry  link  included  a  double-armored , 
0.01-m-diam  electro-mechanical  cable  from  the  transducer  to  a  WES-designed 
Remote  Terminal  Unit  (RTU),  and  a  modem  (Hayes  Smart-modem  1200)  connected 
to  a  commercial  telephone  line.  Data  were  retrieved  daily  through  automatic  dial-up 
and  connection  by  the  mainframe  computer  at  WES.  The  RTU  contained  a  256-Kb 
buffer  to  hold  data  for  short-duration  (3  days)  interruptions  in  the  recovery.  In  the 
event  of  longer-term  failure  of  the  telemetry  link,  a  12-Mb  multi-channel  data  logger 
(Woods  Hole  Instruments)  was  located  onshore.  The  RTU  and  sensors  were  nor¬ 
mally  powered  by  1 10-V  AC  shore  power.  An  independent  battery  pack  would 
operate  the  sensors  and  data  logger  (but  not  the  RTU  and  modem)  for  6  months  in 
the  event  of  power  loss. 

Installation 


The  system  was  installed  during  10  through  13  January  1994.  The  incident 
gauge  (HI01)  was  clamped  about  0.3  m  (1  ft)  off  the  bottom  onto  a  0.05-m-  (0.17- 
ft)-diam  pipe  pile  that  had  been  jetted  3  m  into  the  sand  bottom  in  approximately 
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21m  (68.9  ft)  of  water.1  The  gauge  was  positioned  by  divers  relative  to  the  toe  of 
the  southernmost  extent  of  the  breakwater  head  (range  -  24.4  m  (80  ft);  bearing  - 
30  deg)  and  the  mooring  weight  of  channel  buoy  number  1  (range  -  30.5  m  (100  ft); 
bearing  -  70  deg).  The  interior  gauge  (HI02)  was  clamped  at  the  outer  face  of  the 
wharf  to  pile  No.  23,  or  36.9  m  (121  ft)  east  from  the  west  end  of  the  wharf.  The 
water  depth  in  front  of  the  wharf  was  about  6.4  m  (21  ft),  and  the  transducer  was 
placed  about  1.5  m  (4.9  ft)  above  the  bottom  (Figure  4). 

The  armored  cables  from  the  transducers  were  hand-placed  by  divers  along  the 
toe  of  the  breakwater  and/or  adjacent  to  the  outer  row  of  pilings,  specifically  to 
avoid  disturbance  by  marine  operations.  The  cable  is  sufficiently  dense  that  it  tends 
to  self-scour  several  inches  into  a  sand  bottom  after  a  few  hours.  Both  cables  exited 
the  water  at  the  eastern  end  of  the  wharf.  As  protection  against  topside  operations, 
the  cables  were  encased  in  0.05-m-  (0. 17-ft-)  diameter  schedule  40  galvanized  pipe 
conduit  from  below  the  waterline,  across  the  wharf,  up  the  vertical  rock  cliff  at  the 
east  end  of  the  harbor,  and  across  the  top  of  the  cliff  into  an  existing  compres¬ 
sor/utility  building  (Figure  4).  The  RTU  and  other  electronics  were  mounted  inside 
this  building  in  a  weatherproof  enclosure.  Power  and  telephone  connections  were 
made  at  the  building. 

Results 

Data  collection  commenced  on  16  January  1994,  and  continued  through  Septem¬ 
ber  1994.  The  system  failed  on  29  March  when  a  short  circuit  in  the  power  supply 
to  the  sensors  caused  a  protective  fuse  to  blow.  The  short  occurred  when  the  cable 
to  HI01  was  cut  during  harbor  operations  by  a  barge  mooring  weight.  The  system 
was  diagnosed  and  HI02  brought  back  on  line  on  8  April;  the  cable  to  HI01  was 
repaired  on  28  April.  The  system  operated  essentially  without  interruption  through 
September.  Additional  funds  for  continued  operation  were  unavailable,  so  the 
gauges  were  placed  in  an  inactive  mode  -  data  were  collected,  but  not  analyzed;  no 
data  quality  assurance  was  conducted,  and  no  system  maintenance  was  performed. 
The  system  ceased  operation  on  8  December  1994. 

Summaries  of  significant  wave  height  and  peak  period  are  provided  as  monthly 
plots  in  Appendix  A  (Figures  A1  through  A 18).  Energy  spectra  and  reduced  wave 
parameters  from  each  observation  are  saved  in  a  database  and  are  available  to  POD 
and  WES  researchers  upon  request. 

Wave  conditions  were  not  very  energetic  over  the  course  of  the  study.  The  mean 
incident  wave  height  was  0.5  m  (1.6  ft).  The  largest  wave  occurred  in  January,  with 
a  height  of  1.5  m  (4.9  ft)  and  a  peak  period  of  15.1  sec.  Tables  1  through  4 
summarize  the  data  from  the  two  gauges.  The  deepwater  directional  wave  gauge, 
funded  by  another  project,  was  installed  about  22  km  (13.7  miles)  west  of  Lanai  at 


'Though  mean  depths  can  be  obtained  from  the  pressure  signals,  accurate  depth  data  require 
correction  of  the  measurements  with  atmospheric  pressure  data.  Adjustment  of  the  depth 
information  to  a  common  datum,  such  as  mllw  or  National  Geodetic  Vertical  Datum  involves 
significant  additional  analysis  and/or  measurements.  This  report  does  not  include  water  depth 
information. 
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20.45  deg  N,  157. 13  deg  W  on  6  December  1994,  providing  just  2  days  of  overlap 
before  the  system  ceased  operation.  The  data  logger  for  the  harbor  gauges  has  been 
removed,  but  the  remainder  of  the  system  is  still  in  place. 
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Table  2 

Percent  Occurrence  as  a  Function  of  Peak  Wave  Period  for  Field  Wave 

Gauge  HI01 

Percent  Occurrence  (xIOOO)  of  Height  and  Period 

Height 

m(ft) 

Peak  Period 
secs 

< 

6.9 

6.9- 

8.0 

8.1- 

8.7 

8.8- 

9.5 

9.6- 

10.5 

10.6- 

11.7 

11.8- 

13.3 

13.4- 

15.3 

15.4- 

18.1 

18.2 

< 

Total 

0.0-0.4 
(0.0-1 .4) 

- 

747 

448 

1271 

2767 

6058 

14136 

17352 

3365 

2767 

48911 

0.5-0.9 
(1 .5-3.1) 

- 

224 

373 

1047 

3141 

3964 

13537 

21316 

3365 

1869 

48836 

1.0-1 .4 
(3.2-4.7) 

- 

- 

- 

74 

224 

299 

149 

747 

448 

299 

2240 

1.5-1 .9 
(47-6.4) 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

n 

2.0-2.4 

(6.5-8.0) 

- 

- 

- 

- 

- 

- 

- 

- 

- 

n 

2.5-2.9 
(8.1 -9.6) 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

3.0-3.4 

(9.7-11.3) 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

n 

3.5-3.9 

(11.4-12.9) 

- 

- 

- 

- 

- 

- 

- 

- 

- 

4.0-4.4 

(13.0-14.6) 

- 

- 

- 

- 

- 

- 

- 

- 

- 

n 

4.5-4.9 

(14.7-16.2) 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

0 

5.0+ 

(16.3+) 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

0 

Total 

0 

971 

821 

2392 

6132 

10321 

27822 

39415 

7178 

4935 

Mean  H^,  m  (ft)  =  0.5  (1 .6)  Mean  Tp,  sec  =  13.8 

Highest  H„  m  (ft)  =  1 .5  (4.9)  T otal  Cases  =  1 337 
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Table  3 

Field  Wave  Gauge  HI02  -  Statistical  Summary  of  Wave  Parameters 


Mean  by  Month  - 1994 
m  (ft) 


Mar 

Apr 

May 

0.4 

(1.3) 

0.4 

(1.3) 

0.4 

(1.3) 

Highest  by  Month  - 1994 
m(ft) 


Jan 

Feb 

Mar 

Apr 

May 

Jun 

Jul 

Aug 

Sep 

Oct 

Nov 

Dec 

1.1 

(3.6) 

1.0 

(3.3) 

0.8 

(2.6) 

0.6 

(2.0) 

0.8 

(2.6) 

0.6 

(2.0) 

0.8 

(2.6) 

1.2 

(3.9) 

0.6 

(2.0) 

- 

■ 

- 

One-Year  Statistics  for  Lanai,  HI  (NEMO,  HI01)  (20.79N  156.99W) 


Mean  Hm,  m  (ft) 


Mean  Tp,  sec 


Standard  deviation  of  m  (ft) 


Standard  deviation  of  Tp,  sec 


Highest  m  (ft) 


Tp  associated  with  highest  H^,  sec 


Date  of  highest  occurrence 


0.4  (1.3) 


14.8 


0.1  (0.3) 


1.2  (3.9) 


11.1 


94082420 
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Objectives  and  Approach 


The  numerical  model  study  had  four  main  objectives: 

a  Assist  in  optimizing  the  placement  of  field  wave  gauges. 

b.  Validate  the  model  with  field  data. 

c.  Maximize  efficiency  of  physical  model  experiments  by  identifying  the  most 
promising  harbor  modification  plans  and  most  damaging  incident  wave 
directions. 

d.  Evaluate  harbor  oscillation  characteristics  of  the  final  WES-recommended 
plan  relative  to  the  existing  harbor. 

The  first  objective  was  met  with  a  fast-track  modeling  effort  which  evaluated  the 
existing  harbor  response  to  both  short  waves  (periods  between  6  and  22  sec)  and 
long  waves  (periods  between  23  and  500  sec).  Once  field  gauges  were  in  place  and 
operating,  field  records  with  relatively  high  wave  energy  were  used  for  model 
validation.  Following  validation,  nine  alternative  harbor  plans  were  defined  and 
investigated.  Upon  selection  of  the  final  harbor  plan,  harbor  response  to  long  waves 
was  reevaluated,  including  comparison  with  existing  conditions. 

The  numerical  model,  HARBD,  is  the  standard  WES  tool  for  numerical  harbor 
wave  studies.  The  model  includes  the  following  assumptions: 

a.  No  wave  transmission  through  the  breakwater. 

b.  No  wave  overtopping  of  structures. 

c.  Structure  crest  elevations  above  the  water  surface  cannot  be  studied  or 
optimized. 
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d.  Currents  in  the  channel  cannot  be  evaluated. 

e.  Wave  breaking  effects  in  the  entrance  and  harbor  cannot  be  considered. 

f.  No  nonlinear  effects  are  considered. 

g.  Diffraction  around  structure  ends  is  represented  by  diffraction  around  a  blunt 
vertical  wall  with  specified  reflection  coefficient. 

Despite  limitations  imposed  by  the  above  assumptions,  HARBD  is  considered 
suitable  for  meeting  the  numerical  modeling  objectives  of  the  Kaumalapau  Harbor 
study. 

The  harbor  wave  response  model  is  presented  in  the  following  section,  including 
a  general  description  of  the  HARBD  model  and  implementation  of  the  model  at 
Kaumalapau  Harbor.  Validation  was  accomplished  with  eight  high  wave  cases 
selected  from  available  field  data.  The  final  section  of  this  chapter  describes  the 
numerical  simulation  procedures  and  calculations. 

As  part  of  the  simulation  procedures,  a  suite  of  incident  wave  conditions  must  be 
specified  at  the  seaward  boundary  of  the  area  covered  by  HARBD.  Often  the 
incident  short  waves  are  determined  by  consideration  of  offshore  wave  climate  and 
wave  transformation  when  propagating  over  nearshore  bathymetry.  This  step  was 
not  included  in  the  Kaumalapau  Harbor  study.  Instead,  the  full  range  of  possible 
incident  wind  wave  and  swell  conditions  was  represented  in  the  incident  waves 
studied.  A  broad  range  of  long  wave  conditions  was  also  studied  to  identify 
possible  harbor  resonant  responses. 

Results  for  wind  waves  and  swell  are  presented  in  Chapter  5.  Harbor  oscillation 
results  are  presented  in  Chapter  6.  The  presentation  focuses  on  wave  conditions  in 
the  vicinity  of  the  main  dock,  but  results  over  the  full  harbor  area  also  are  given. 


Model  Description 


Model  formulation 

The  numerical  wave  model  HARBD  is  a  steady-state  hybrid  element  model 
used  in  the  calculation  of  linear  wave  response  in  harbors  of  varying  size  and 
depth  (Chen  1986,  Chen  and  Houston  1987,  Lillycrop  and  Thompson  1995). 
Originally  developed  for  use  with  long-period  waves  (Chen  and  Mei  1974), 
HARBD  has  since  been  adapted  to  include  capabilities  for  modeling  wind  waves 
and  swell  (Houston  1981),  bottom  friction,  and  partially  reflective  boundaries 
(Chen  1986).  The  model  is  based  on  a  linearized  mild  slope  equation.  An 
overview  of  the  model  and  its  applications  is  given  by  Thompson  and  Hadley 
(1995). 

The  HARBD  model  has  been  shown  to  perform  satisfactorily  in  comparison  to 
analytic  solutions  and  laboratory  data  for  a  variety  of  wind  wave  and  swell  cases 
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(Houston  1981;  Crawford  and  Chen  1988;  Thompson,  Chen,  and  Hadley  1996) 
and  long-wave  cases  (Chen  1986;  Chen  and  Houston  1987;  Houston  1981; 
Thompson,  Chen,  and  Hadley  1993).  As  a  result,  it  has  been  used  with  confi¬ 
dence  in  both  long-wave  and  short-wave  studies.  Long-wave  studies  have 
included  harbor  oscillations  (Thompson  and  Hadley  1994b;  Briggs  et  al.  1994; 
Briggs,  Lillycrop,  and  McGehee  1992;  Mesa  1992;  Sargent  1989;  Weishar  and 
Aubrey  1986;  Houston  1976)  and  tsunamis  (Farrar  and  Houston  1982,  Houston 
and  Garcia  1978,  Houston  1978).  Wind  wave  and  swell  studies  include  Thomp¬ 
son  and  Hadley  (1994a,  1994b);  Lillycrop  et  al.  (1993);  Lillycrop  and  Boc 
(1992);  Lillycrop,  Bratos,  and  Thompson  (1990);  Kaihatu,  Lillycrop,  and 
Thompson  (1989);  Farrar  and  Chen  (1987);  Clausner  and  Abel  (1986);  and 
Bottin,  Sargent,  and  Mize  (1985). 

The  HARBD  model  covers  in  detail  a  domain  including  the  harbor  and  a 
portion  of  the  adjacent  nearshore  area  (Figure  5).  This  domain  is  bounded  by  a 
180-deg  semicircle  in  the  water  region  seaward  of  the  harbor  entrance  (d A  in 
Figure  5)  and  the  land-water  interface  along  the  shoreline  and  harbor  (5C  in 
Figure  5).  The  region  defined  by  these  boundaries  is  denoted  Region  A.  If 
possible,  the  semicircle  radius  should  be  at  least  twice  the  wavelength  of  the 
longest  incident  wave  to  be  modeled  (using  a  typical  water  depth  within  the 
semicircle).  Also,  the  semicircle  should  encompass  any  complex  offshore 
bathymetry  which  strongly  influences  waves  entering  the  harbor.  In  general,  the 
semicircle  should  be  as  large  as  practical  constraints  on  grid  size  and  resolution 
will  allow. 

The  area  outside  the  semicircle  is  treated  as  a  semi-infinite  region  which 
extends  from  a  straight  coastline  seaward  to  infinity  (Region  B).  This  region  is 
assumed  to  have  a  constant  water  depth  and  no  bottom  friction. 

Assuming  linear,  regular  waves  propagating  over  mild  slope  in  arbitrary  water 
depth,  Chen  (1986)  derived  the  governing  equation  as 

V  •  (  X  c  c  V<f>  )  +  OLft  c()  =  0  (1) 

s  c 


where 

V  =  horizontal  gradient  operator1 
A  =  complex  bottom  friction  factor 
c  =  wave  phase  speed 
Cj  =  wave  group  speed 
<p  —  velocity  potential 
co  —  angular  frequency 

This  equation  is  identical  to  Berkhoff  s  (1972)  equation  except  for  addition  of  the 
bottom  friction  factor  A.  The  factor  A,  which  is  a  complex  number  with  magni¬ 
tude  greater  than  zero  and  less  than  or  equal  to  one,  is  specified  as 


i 


For  convenience,  symbols  and  abbreviations  are  listed  in  the  notation  (Appendix  C). 
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Figure  5.  Representation  of  HARBD  domain 


X 


d  sinh  k d 


j'Y 


(2) 


where 

i  =(-l)1/2 

P  =  dimensionless  bottom  friction  coefficient  that  can  vary  in  space 
a,  =  incident  wave  amplitude 
d  =  water  depth 
k  =  wave  number 

Y  =  phase  shift  between  stress  and  flow  velocity 


The  bottom  friction  factor  is  a  factor  tending  to  reduce  local  velocities  propor¬ 
tionately  through  the  relationships 


u 


V 


(3) 
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where 


u,v  =  local  horizontal  velocity  components 
x ,y  =  horizontal  coordinates 

Boundary  conditions  are  specified  in  Regions  A  and  B.  At  the  solid  boundary 
dC,  a  reflection/absorption  boundary  condition  is  used  similar  to  the  impedance 
condition  in  acoustics.  The  condition  is  specified  as 

fr -  “  *  *  0  M 


with 


a 


ZK 


1  -K 
1  +  Kr 


(5) 


where 

n  =  unit  normal  vector  directed  into  the  solid  region 

Kr  —  reflection  coefficient  of  the  boundary 

Values  of  Kr  for  wind  waves  and  swell  are  normally  chosen  based  on  the 
boundary  material  and  shape.  General  guidelines  for  Kr  can  be  assembled  from 
laboratory  and  field  data  (Thompson,  Chen,  and  Hadley  1996).  In  wind  wave  and 
swell  studies,  Kr  is  generally  chosen  to  be  consistent  with  this  guidance.  Effects 
such  as  slope,  permeability,  relative  depth,  wave  period,  breaking,  and  overtop¬ 
ping  can  be  considered  in  selecting  values  within  these  fairly  wide  ranges.  For 
long-wave  studies,  Kr  is  generally  set  equal  to  1.0,  representing  full  reflection. 

The  second  boundary  condition  is  imposed  in  the  far  region  (Region  B )  at 
infinity.  It  requires  that  the  scattered  wave,  defined  as  the  difference  between  the 
total  wave  and  incident  wave,  behave  as  a  classical  outgoing  wave  at  infinity. 

This  radiation  condition  may  be  expressed  as 

Um  ^  -  *)  <0-  -  0  (6) 

r-  °° 


where 

r  =  radial  polar  coordinate 

4?  =  velocity  potential  of  the  scattered  wave 

d  =  partial  differential  symbol 


The  complete  boundary  value  problem  is  specified  by  Equations  1,  4,  and  6. 
A  hybrid  element  method  is  employed  to  solve  the  boundary  value  problem.  A 
conventional  finite  element  grid  is  developed  and  solved  in  Region  A.  The 
triangular  elements  allow  detailed  representation  of  harbor  features  and 
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bathymetry  within  Region  A.  An  analytical  solution  with  unknown  coefficients  in 
a  Hankel  function  series  is  used  to  describe  Region  B.  For  a  given  grid,  short 
wave  period  simulations  (relatively  large  values  of  at)  require  more  terms  than 
long-period  simulations  to  adequately  represent  the  series.  A  variational  principle 
with  a  proper  functional  is  established  such  that  matching  conditions  are  satisfied 
along  <34.  Details  are  given  by  Chen  (1986)  and  Lillycrop  and  Thompson  (1995). 

Experience  with  the  model  has  indicated  that  the  element  size  Ax  and  local 
wavelength  L  should  be  related  by 

A*  *  §  (7) 


Typically,  harbor  domains  include  some  shallow  areas  in  which  many  elements 
would  be  needed  to  satisfy  the  constraint  in  Equation  7.  In  practice,  Equation  7  is 
at  least  satisfied  in  the  harbor  channel  and  basin  depths.  If  additional  elements 
can  be  accomodated,  it  is  generally  preferred  to  extend  the  semicircle  further 
seaward  rather  than  to  greatly  refine  shallow  harbor  regions. 

Input  information  for  HARBD  must  be  carefully  assembled.  In  addition  to 
developing  the  finite  element  grid  to  suit  HARBD  requirements,  a  number  of 
parameters  must  be  specified.  Critical  input  parameters  and  ranges  of  typical 
values  are  summarized  in  Table  5. 


Table  5 

Critical  HARBD  Input  Parameters  and  Ranges  of  Typical  Values 

Typical  Values 

Parameter 

Where  Specified 

Short  Waves 

Long  Waves 

Bottom  friction,  P 

Every  element 

0.0 

0.0  -  0.1 

Boundary  reflection,  K, 

Every  element  on  solid  boundary 

0.0  -  1 .0 

1.0 

Coastline  reflection,  Krcm , 

Single  value 

1.0 

1.0 

Depth  in  infinite  region, 

diet 

Single  value 

Between  avg.  &  max.  on  semicircle 

Number  of  terms  in 

Hankel  function  series 

Single  value 

8  -  100’ 

8 

’The  number  of  terms  needed  increases  as  wave  period  decreases. 

The  principal  output  information  available  from  HARBD  consists  of  amplifica¬ 
tion  factor  and  phase  at  each  node.  These  are  defined  as 
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(8) 


amp 


!4>l 


0 


tan"1 


Im  {())} 
Re  {<|>} 


where 

=  amplification  factor 

a,  Of  =  local  and  incident  wave  amplitudes 
H,  Hi  =  local  and  incident  wave  heights 
6  =  phase  relative  to  the  incident  wave 
Im{<p}  =  imaginary  part  of  (p 
Re{ <p}  =  real  part  of  <p 


Amplification  factors  are  easily  interpreted.  Phases  are  helpful  in  viewing  wind 
wave  and  swell  propagation  characteristics  and  in  interpreting  standing  wave 
patterns.  In  long-wave  applications,  phases  prove  useful  for  determining  relative 
phase  differences  within  the  harbor,  interpreting  harbor  oscillation  patterns,  and 
identifying  potentially  troublesome  nodal  areas. 

Spectral  adaptation 


Often  the  model  is  used  to  approximate  irregular  wind  wave  and  swell  behavior, 
as  in  physical  model  experiments  with  irregular  waves  and  all  field  cases.  More 
realistic  numerical  model  simulations  can  be  obtained  by  linearly  combining 
HARBD  results  from  a  range  of  regular  wave  frequencies  in  the  irregular  wave 
spectrum.  Regular  wave  results  are  weighted  to  properly  represent  the  desired 
spectral  distribution  of  energy.  The  concept  of  linear  superposition  of  weighted 
regular  wave  results  also  can  be  extended  to  include  directional  spreading  in  the 
spectrum  to  be  simulated. 

Spectral  adaptation  of  the  HARBD  model  is  done  as  a  post-processing  step 
using  the  standard,  regular  wave  output  from  the  model.  For  a  given  incident  wave 
direction,  HARBD  is  run  for  a  number  of  wave  periods  spread  between  the  shortest 
period  satisfying  the  grid  resolution  constraint  of  Equation  7  and  the  longest  swell 
period  of  interest. 


Spectral  post-processing  is  based  on  the  assumption  that  a  consistent  spectral 
form  can  be  applied  at  every  node.  This  major  assumption  provides  the  basis  for  a 
workable,  reasonable  spectral  weighting  which  improves  on  the  traditional  regular 
wave  approach.  The  JONSWAP  spectral  form  was  chosen  (Hasselmann  et  al. 
1973).  The  JONSWAP  spectrum  is  specified  as  (U.S.  Army  Corps  of  Engineers 
1989) 


S(f) 


«  g2  wb 

(2*  ff- 


(9) 
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where 


S(fi)  =  spectral  energy  density  at  frequency  f 

The  parameters  a  and  b  are  given  by  the  following  relationships 

-1.25 


a  = 


f  Tn 

j  i  p 


b  =  e 


2<J' 


a  =  0.07  for  /.</ 


=  0.09  for  fef 


where 


Tp  =  peak  spectral  period  j 

fp  =  peak  spectral  frequency  =  — 


Parameters  a:  and  y  are  calculated  as 

a  =  157.9  e2 

y  =  6614  e159 

e  .  h- 

4  L 

p 


(10) 


(11) 


where 

Hs  =  significant  wave  height 

Lp  =  wavelength  for  waves  at  peak  frequency 

The  parameter  eis  a  significant  wave  steepness.  The  parameter  y,  called  the  peak 
enhancement  factor,  controls  the  sharpness  of  the  spectral  peak. 
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Although  the  JONSWAP  spectrum  was 
developed  primarily  for  actively  growing 
wind  waves,  it  can  be  used  to  approximate 
any  single-peaked  spectrum,  including  old 
swell  which  has  travelled  a  great  distance 
from  the  generation  area,  by  appropriate 
choice  of  y  (e.g.  Goda  1985)  (Table  6). 

Spectral  post-processing  begins  with 
specification  of  the  desired  Hp  Tp  and  y 
and  the  arrays  of  HARBD  amplification 
factors.  A  refined  JONSWAP  spectrum  is 
computed  with  1,000  points,  where  the  f/s  in  Equation  9  are 

f  =  0.5% ,  f2  =  0.502 % ,  f3  =  0.504% , ... ,  f1000  =  2.498% 

The  number  of  wave  periods  computed  with  HARBD  is  always  much  smaller 
than  1,000,  typically  less  than  20.  These  periods,  converted  to  frequency  (reciprocal 
of  period),  can  be  used  to  define  bands  in  the  JONSWAP  spectrum.  Bands  are 
bounded  by  the  midpoints  between  HARBD  computational  frequencies.  The 
highest  and  lowest  frequency  bands  are  assumed  to  be  centered  on  the  highest  and 
lowest  HARBD  computational  frequencies,  respectively.  A  weighting  factor  for 
each  HARBD-defmed  band  is  computed  by  summing  values  from  the  refined 
JONSWAP  spectrum  which  fall  within  the  band  and  normalizing  by  the  total 
spectral  energy. 

E  m 

i'Ntl 

=  "iooo - 

E  S(f) 

i=  1 


(12) 


where 

wk  =  weighting  factor  for  &’th  HARBD  computational  frequency 

Nu  =  index  of  highest  JONSWAP  frequency,/;,  satisfying  f.  < 

fk-ifbfk+i  =  (fc-l)’th,  £’th,  and  (&+l)’th  HARBD  computational  frequencies, 
with  fk.,<fk<fk+1 

Though  not  shown  in  the  equation,  the  weighting  factor  also  includes  fractional 
energy  interpolated  across  JONSWAP  frequencies  bracketing  the  two  end  points  of 
each  HARBD  band. 
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The  effective  amplification  factor  at  each  node  is  computed  as 


t 


amp'eff 


(13) 


where 

(AamP),ff  -  effective ,  or  spectral,  amplification  factor  at  a  node 
A  (fj  =  nodal  amplification  factor  for  HARBD  computational  frequency  y* 
NT  =  number  of  HARBD  computational  wave  periods 


Finite  element  grids 

The  finite  element  numerical  grid  depicting  existing  conditions  at  Kaumalapau 
Harbor  was  created  using  WES's  finite  element  grid  development  software  (Turner 
and  Baptista  1993)  (Figure  6).  The  grid  covers  the  entire  Kaumalapau  Harbor 
embayment  and  extends  somewhat  seaward  from  the  bay  entrance.  The  land 
boundary  was  digitized  from  a  NOAA  nautical  chart.  Grid  element  size  is  based  on 
the  criteria  of  six  elements  per  wavelength  (the  minimum  recommended  resolution 
with  HARBD)  for  a  6-sec  wave  in  a  2.7-m  (9-ft)  water  depth.  Depths  over  virtually 
the  entire  embayment  exceed  2.7  m  (9  ft).  For  the  longer  period  waves,  the  grids 
give  a  high  degree  of  resolution.  Grid  characteristics  are  summarized  in 
Table  7. 

The  radius  of  the  seaward  semicircle  is  303  m  (995  ft).  This  is  equivalent  to  0.9 
and  5.4  wavelengths  for  the  longest  and  shortest  short  wave  periods  considered,  as¬ 
suming  a  representative  water  depth  of  27.4  m  (90  ft).  The  semicircle  size  and 
location  were  chosen  to  include  the  breakwater,  the  large  shoal  north  of  the  break¬ 
water,  and  the  point  and  shoal  area  south  of  the  harbor  entrance.  The  semicircle 
extends  sufficiently  far  seaward  to  cover  the  most  important  nearshore 
bathymetry. 

Bathymetric  data  were  obtained  from  NOAA  hydrographic  chart  19351. 
Digitized  depths  were  transferred  onto  the  finite  element  grid  using  the  WES  grid 
software  package.  A  contour  plot  of  bathymetry  is  given  in  Figure  7. 

Reflection  coefficients  Kr  are  needed  for  all  solid  boundaries.  For  the  short  wave 
simulations,  Kr  values  were  estimated  from  existing  Corps  of  Engineers  guidance, 
photos  from  a  recent  site  visit  by  WES  personnel,  and  past  experience.  The  solid 
boundary  was  divided  into  eight  zones  and  a  reflection  coefficient  was 
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Figure  6.  Grid  of  existing  harbor 


Table  7 

Grid  Sizes 

Breakwater  Layout 

Number  of: 

Length 
of  Typical 
Element 
m(ft)  I 

Elements 

Nodes 

Solid 

Boundary 

Nodes 

Semicircle 

Boundary 

Nodes 

Existing 

17,376 

8,929 

285 

196 

4.9  (16.1) 

Straight  Extension 

17,364 

8,927 

296 

196 

4.9  (16.1) 

Dogleg  Extension 

17,355 

8,924 

296 

196 

4.9  (16.1) 

Southern  Breakwater 

17,350 

8,927 

307 

196 

4.9  (16.1) 

WES-recommended 

Plan 

17,396 

8,945 

297 

196 

4.9  (16.1) 
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Depths  in  Meters 


Figure  7.  Bathymetry 


estimated  for  each  zone  (Figure  8).  Reflection  coefficients  ranged  from  0.2  for  the 
shallow  sandy  area  around  the  stream  outlet  to  0.5  for  the  exposed  cliff  areas  along 
much  of  the  embayment  to  0.8  along  the  wharf  face.  Other  parameter  values  used  in 
the  numerical  model  are  summarized  in  Table  8. 

Different  parameters  are  used  for  the  long-wave  simulations.  The  reflection  co¬ 
efficient  was  set  to  1.0  for  all  boundaries,  since  long  waves  generally  reflect  very 
well  from  a  coastal  boundary.  Long  waves  are  more  affected  by  bottom  friction 
than  short  waves,  so  a  value  of  /?  greater  than  zero  is  appropriate.  The  value  of  p  is 
best  determined  by  calibration  with  field  data.  However,  data  records  in  this  study 
are  too  short  to  provide  reliable  long-wave  information.  An  accurate  value  for  /?  is 
not  critical  to  objectives  of  the  study,  and  a  default  value  of  /?=  0  was  used.  This 
and  other  parameters  are  summarized  in  Table  8. 
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Figure  8.  Wave  reflection  coefficient  values,  short  waves 


Table  8 

Parameter  Values  Used  in  HARBD 

Parameter 

Value 

Short  Waves 

Long  Waves 

Bottom  friction,  /? 

0.0 

0.0 

Coastline  reflection,  Krr^« 

1.0 

1.0 

Depth  in  infinite  region,  d,„ 

38.4  m  (126  ft) 

38.4  m  (126  ft) 

Nine  plans  for  modifying  the  harbor  were  defined  after  a  study  team  review 
meeting  with  sponsors  (Table  9).  The  plans  involved  three  breakwater  modifica¬ 
tions,  including  a  breakwater  at  the  southern  end  of  the  harbor  embayment  and  two 
alternatives  for  extending  the  existing  breakwater.  During  physical  model  experi¬ 
ments,  a  final  WES-recommended  plan  was  developed.  The  existing  harbor  grid 
was  modified  to  represent  the  alternative  breakwater  configurations  (Figures  9 
through  12).  Grid  characteristics  for  each  of  the  four  breakwater  configurations  are 
included  in  Table  7. 
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Table  9 

Harbor  Alternatives  for  Numerical  Modeling 


Feature 

Test  Series 

North  breakwater 

Existing 

Extended  (straight) 

Extended  (dogleg) 

X 

X 

X 

X 

X 

X 

X 

1 

X 

South  breakwater 

X 

■ 

Wave  absorber 

Northeast  shore  (full  length) 

Northeast  shore  (north  half  only) 

South  shore  (full  length) 

X 

X 

X 

X 

1 

1 

Figure  9.  Grid  with  dogleg  breakwater  extension 
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Figure  10.  Grid  with  straight  breakwater  extension 


Figure  1 1 .  Grid  with  south  breakwater  addition 
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Figure  12.  Grid  of  WES-recommended  plan 


The  harbor  modifications  studied  include  combinations  of  the  following: 

a.  45.7-m  (150-  ft)  dogleg  extension  to  the  existing  breakwater,  angled  into  the 
harbor. 

b.  45 . 7-m  ( 1 5 0-ft)  straight  extension  to  the  existing  breakwater. 

c.  6 1  -m  (200-ft)  rubble-mound  breakwater,  extending  northwest  from  the 
southern  comer  of  the  harbor  entrance. 

d.  Addition  of  wave  absorber  along  the  northeast  portion  of  the  harbor  shore. 

e.  Addition  of  wave  absorber  along  the  northernmost  portion  of  the  northeast 
harbor  shore. 

f.  Addition  of  wave  absorber  along  the  southeast  portion  of  the  harbor  shore. 

The  addition  of  wave  absorber  was  simulated  by  reducing  the  value  of  the  short¬ 
wave  reflection  coefficient  to  Kr  =  0.2  over  the  affected  area.  Absorber  was  simu¬ 
lated  over  large  segments  of  the  harbor  boundaiy  to  assess  the  maximum  potential 
effect  of  this  type  of  modification  (Figure  13). 
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Figure  1 3.  Wave  reflection  coefficient  values  for  wave  absorber 
simulations,  short  waves 


Validation 

The  availability  of  field  data  at  Kaumalapau  Harbor  allowed  a  validation  test  of 
the  numerical  harbor  response  model.  Validation  was  confined  to  short-wave  cases. 
The  field  records  were  not  suited  to  long-wave  validation,  and  validation  was  not 
essential  for  meeting  the  comparative  objectives  of  the  long-wave  study. 

Field  validation  of  a  short-wave  model  is  inherently  difficult,  typically  including 
uncertainties  about  some  aspects  of  incident  waves,  bathymetry,  and  construction 
underneath  wharves  and  piers  and  inside  breakwaters  (which  affects  harbor  bound¬ 
ary  locations  and  reflection  coefficients).  Even  when  the  structure  is  known,  the 
effective  reflection  coefficient  can  rarely  be  estimated  with  precision.  Despite  the 
impossibility  of  getting  a  highly  accurate  and  controlled  field  validation  experiment, 
field  data  are  valuable  for  a  rough  validation  to  ensure  that  the  model,  with  parame¬ 
ter  values  within  the  range  of  reasonable  choices  given  the  particular  site,  gives  re¬ 
sults  consistent  with  measurements. 

Five  events  occurred  during  the  winter  of  1994  in  which  Hm0  at  the  outside  gauge 
exceeded  1  m.  One  or  two  records  within  each  of  these  events  was  selected  for 
HARBD  simulation  (Table  10).  The  main  criteria  for  selection  are  that  Hm0  and  Tp 
at  Kaumalapau  Harbor  and  mean  wave  direction  from  the  NOAA  buoy  north  of  the 
island  of  Molokai  (Gauge  51026)  are  relatively  stable  for  at  least  a  few  hours 
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Table  10 

Summary  of  Field  Cases  for  Model  Validation 

Date 

Hour 

GMT 

Wave  Period 

Sec 

Wave  Direction 

Deq 

17  Jan  94 

1200 

15.1 

300 

31  Jan  94 

0800 

19.7 

325 

31  Jan  94 

1200 

18.3 

325 

24  Feb  94 

1200 

11.1 

320 

27  Feb  94 

0400 

17.1 

301 

27  Feb  94 

0800 

17.1 

301 

11  Mar  94 

1200 

15.1 

300 

12  Mar  94 

1200 

12.8 
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preceding  the  selected  record  and  the  spectral  shape  is  single  peaked  if  possible. 
These  criteria  help  to  ensure  that  Hm0  and  Tp  are  well-defined  and  that  the 
JONSWAP  spectral  form  imposed  on  the  numerical  model  results  is  at  least  qualita¬ 
tively  correct.  All  five  of  the  events  also  appeared  as  high  wave  events  in  the 
Molokai  buoy  measurements.  This  observation,  coupled  with  the  climatological 
predominance  of  winter  storm  events  from  northerly  directions,  strongly  implies  that 
the  waves  at  Kaumalapau  were  coming  from  the  northwest.  Some  exploratory 
HARBD  runs  were  done  using  incident  wave  directions  on  either  end  of  the  reason¬ 
able  bounds  at  Kaumalapau,  consistent  with  a  northerly  approach  and  sheltering  by 
the  northwest  tip  of  Lanai.  Since  results  were  only  mildly  sensitive  to  incident  di¬ 
rection,  the  actual  direction  reported  from  the  Molokai  buoy  (which  fell  within  the 
reasonable  bounds  at  Lanai)  was  used  in  the  final  HARBD  runs.  Frequency-de¬ 
pendent  amplification  factors  from  the  HARBD  runs  were  weighted  based  on  the 
energy  levels  in  a  JONSWAP  spectrum.  Appropriate  Tp  and  spectral  shape  parame¬ 
ters  were  used  to  establish  the  weightings. 

Reflection  coefficient  at  the  wharf  is  a  key  model  parameter  in  comparing  to 
measurements  at  mid-wharf.  The  wharf  face  is  a  concrete  deck  supported  by  con¬ 
crete  piles.  Under  the  wharf,  within  about  2  to  3  m  (6.6  to  9.8  ft)  landward  of  the 
wharf  face,  is  a  rough  vertical  wall  of  solid  rock.  It  was  estimated  that  a  reflection 
coefficient  between  0.6  and  0.9  would  be  representative  of  the  wharf.  Initially  Kr  = 
0.9  was  used  to  be  conservative. 

The  ratio  of  Hm0's  at  the  inner  and  outer  harbor  gauges  was  compared  to  amplifi¬ 
cation  factor  ratio  at  the  same  locations  from  HARBD  runs  using  the  spectral 
approach.  The  numerical  model  results  showed  a  tendency  to  overestimate  wave 
conditions  at  the  wharf.  Runs  were  repeated  with  wharf  reflection  coefficient 
adjusted  to  0.8  and  0.6.  A  reflection  coefficient  of  0.8  gave  satisfactory,  unbiased 
results  (Figure  14).  The  figure  also  indicates  the  location  of  the  data  cloud  when  the 
reflection  coefficient  at  the  wharf  is  set  to  0.6  and  0.9,  excluding  for  clarity  the 
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Figure  14.  Model  validation 


single  point  for  which  the  field  ratio  was  greater  than  1.0.  The  best-fitting  reflection 
coefficient,  0.8,  is  very  reasonable  for  this  wharf  structure,  and  the  field  validation  is 
considered  successful. 

Varying  reflection  coefficient  along  the  wharf  significantly  affected  wave  esti¬ 
mates  at  the  wharf,  but  the  effect  was  very  localized.  Areas  away  from  the  immedi¬ 
ate  wharf  area  were  relatively  unchanged.  Similar  localized  changes  appeared  when 
boundary  reflection  coefficients  were  varied  as  part  of  the  short  wave  simulations  of 
harbor  modification  alternatives  (Chapter  5). 


Simulation  Procedures  and  Calculations 

Incident  wave  conditions 

A  range  of  short-and  long-wave  conditions  incident  to  the  Kaumalapau  Harbor 
breakwater  and  embayment  was  considered.  A  full  range  of  wave  periods  and  direc¬ 
tions  which  could  cause  damaging  waves  inside  the  harbor  was  included. 
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The  short  wave  periods  and  ap¬ 
proach  directions  considered  are  given 
in  Table  1 1 .  The  shortest  wave  period 
represents  a  local  storm  condition 
which,  based  on  local  experience,  can 
generate  waves  of  concern  within  the 
harbor.  The  longest  period  represents 
a  very  long  sw'ell  condition.  Directions 
were  chosen  to  include  all  likely  ap¬ 
proach  directions  to  the  harbor  en¬ 
trance.  Numerical  simulation  direc¬ 
tions  were  reckoned  in  10-deg  incre¬ 
ments  relative  to  the  straight-in  ap¬ 
proach  direction  (261.5-deg  azimuth, 
rounded  to  261  deg  for  determining 
calculation  directions  ).  Incident  wave 
directions  and  the  angular  orientation 
of  the  seaward  semicircular  model 
boundary  are  illustrated  in  Figure  15. 

Not  all  of  the  incident  short  wave 
conditions  listed  in  the  table  were  con¬ 
sidered  in  every  set  of  simulations. 

Directional  sensitivity  tests  of  the  ex¬ 
isting  harbor  included  HARBD  calcu¬ 
lations  with  all  of  the  listed  wave  peri¬ 
ods  and  all  of  the  listed  directions. 

Directional  sensitivity  then  was  evalu¬ 
ated  for  spectra  with  Tp  values  of  8,  10, 

16,  and  18  sec,  representative  of  storm  and  swell  conditions  at  the  site. 

For  comparing  alternatives,  HARBD  was  run  with  the  full  set  of  short  wave  peri¬ 
ods  for  each  simulation  and  a  representative  set  of  directions,  chosen  after  consider¬ 
ation  of  the  directional  sensitivity  test  results.  The  numerical  simulation  directions 
were  the  20 1-,  22 1-,  25 1-,  29 1-,  and  3 15 -deg  azimuth.  The  3 15-deg  direction  was 
used  in  place  of  321  deg  to  represent  a  more  realistic  northerly  limit  on  incident 
wave  direction.  Alternatives  were  evaluated  for  spectra  with  Tp  values  of  10  and 
16  sec. 

Wave  heights  and  probabilities  of  occurrence,  which  also  generally  are  needed  in 
harbor  studies,  were  not  necessary  in  this  study.  The  objectives  of  the  study  were  to 
assess  harbor  characteristics  for  various  harbor  layouts,  but  did  not  include  the  esti¬ 
mation  of  wave  climate  inside  the  harbor.  Hence,  an  estimate  of  local  wave  climate 
was  not  essential.  Since  the  numerical  model  used  is  linear,  the  results  for  each  inci¬ 
dent  wave  period  and  direction  can  be  taken  as  representative  for  any  nonbreaking 
wave  height. 
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Figure  15.  Incident  wave  directions 


Incident  long-wave  conditions  considered  are 
given  in  Table  12.  A  fine  resolution  in  wave  fre¬ 
quency  was  used  over  the  full  range  of  possible 
resonant  conditions  to  ensure  that  all  important 
peaks  were  identified.  A  total  of  202  periods  were 
considered.  Only  one  approach  direction  is  in¬ 
cluded,  since  past  studies  have  indicated  that  har¬ 
bor  response  is  relatively  insensitive  to  incident 
long-wave  direction.  This  direction  represents  a 
wave  directly  approaching  the  harbor  entrance 
from  deep  water. 

One  water  level  was  used  in  the  simulations. 
The  tide  range  at  Kaumalapau  Harbor  is  relatively 
small,  with  a  mean  range  of  0.46  m  (1.5  ft).  A 
high-water  condition  corresponding  to  the  mean 
high-water  tide  level  was  selected  because  high- 
wave  conditions  are  often  accompanied  by  some- 


Table  12 

Summary  of 

Incident  Long  Wave 
Conditions 

Wave 

Period 

sec 

Wave  Direction 
deg  azimuth 

22.0 

261.5 

22.3 

22.6 

1 

513.0 

1  Frequency  increments  are 
0.0001  Hz  for  periods  of  22- 
82  sec  and  0.00006  Hz  for 
periods  of  83-513  sec. 
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what  elevated  water  levels.  Mean  high  water  at  Kaumalapau  Harbor  is  +0.52  m 
(+1.7  ft)  mllw. 

Calculation  of  spectra 

All  numerical  model  simulation  results  for  Kaumalapau  Harbor  are  based  on 
spectral  post-processing  of  the  initial  HARBD  runs.  This  approach  requires,  first, 
that  HARBD  be  run  with  the  range  of  wave  periods  to  be  considered  in  the  spectral 
calculations.  Second,  values  of  peak  wave  period  Tp,  corresponding  to  the  peak 
spectral  frequency,  and  spectral  peak  enhancement  factor  y  must  be  specified.  The 
Tp  values  were  chosen  to  represent  wind  wave  and  swell  conditions  at  the  harbor,  as 
discussed  in  the  section  "Incident  Wave  Conditions." 

Values  for  y  were  approximated  by  re¬ 
lating  the  guidance  in  Table  6  to  Tp  values 
(Table  13).  High-energy  waves,  of  concern 
for  harbor  design,  with  Tp  up  to  10  sec  were 
assumed  to  be  growing  seas,  and  y  was  set 
to  3.3.  Waves  with  Tp  equal  to  16  sec  and 
18  sec  were  treated  as  swell.  As  swell  Tp 
increases,  the  swell  is  expected  to  have  an 
increasingly  peaked  spectrum.  However,  a 
value  of  y  =  8  was  considered  representa¬ 
tive  of  the  two  swell  periods  considered. 


Output  basins 

In  order  to  get  special  coverage  of  areas 
where  harbor  traffic  would  most  likely  be 
affected  by  wave  conditions,  1 1  output  loca¬ 
tions  or  "basins"  were  selected  in  the  ap¬ 
proach  to  the  main  wharf  (basins  5  through  7  and  1 1),  the  main  wharf 
(basins  lthrough  4),  and  peripheral  areas  (basins  8  through  10)  (Figure  16).  Basins 
1  and  5  coincide  with  field  wave  gauge  locations.  A  basin  is  a  small  cluster  of  ele¬ 
ments  over  which  the  HARBD  response  is  averaged  to  give  a  more  representative 
output.  Each  basin  in  this  study  contains  41  to  54  elements.  At  each  of  these  loca¬ 
tions,  HARBD  output  information  was  saved  in  addition  to  the  detailed  output  at 
nodes. 


Table  13 

Approximate 

Relationship 

Between  T„  and  y 

TP 

Y 

sec 

8 

3.3 

10 

3.3 

16 

8 

18 

8 
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5  Harbor  Response  to  Wind 
Waves  and  Swell1 


Numerical  model  studies  of  the  harbor  response  to  wind  waves  and  swell  were 
directed  toward  increasing  the  efficiency  and  decreasing  the  time  and  cost  for 
physical  model  experiments. 

Directional  Sensitivity  of  the  Harbor 

Large  wavemakers  were  used  to  generate  unidirectional  irregular  waves  in  the 
physical  model.  To  change  the  incident  wave  direction,  the  wavemakers  had  to 
be  physically  moved  and  side  training  walls  relocated,  a  time-consuming  and 
expensive  operation.  For  practical  reasons,  only  three  incident  directions  could 
be  considered.  To  assure  that  the  directions  were  representative  of  the  most 
damaging  harbor  conditions,  a  directional  sensitivity  study  was  performed  with  the 
HARBD  model  for  incident  Tp  values  of  8,  10,  16,  and  18  sec  (Figures  17  through 
20). 

Figures  17  and  18  indicate  that  shorter  period  waves  at  the  wharf  have  a 
tendency  to  be  high  (relative  to  incident  waves)  when  coming  from  a  southerly 
direction,  though  the  directional  response  for  8-sec  waves  also  peaks  around  a 
northerly  approach  direction  of  281  deg.  Since  the  wharf  is  most  directly  exposed 
to  waves  from  the  southwest,  a  strong  response  to  southerly  waves  is  not  surpris¬ 
ing.  For  the  longer  period  swell,  the  peak  directional  response  shifts  more  toward 
the  direct  approach  to  the  harbor  of  251  to  261  deg  (Figures  19  and  20). 

Comparison  of  Alternatives 

Nine  different  harbor  modification  alternatives  were  defined  by  the 
Kaumalapau  Harbor  study  team  for  exploratory  investigation,  as  discussed  earlier. 
Alternatives  were  assessed  with  the  numerical  model.  Physical  model  experiment 
plans  were  streamlined  based  on  conclusions  from  the  numerical  model  simula¬ 
tions. 
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Figure  1 8.  Amplification  factor  sensitivity  to  wave  direction,  Tp 
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Figure  19.  Amplification  factor  sensitivity  to  wave  direction,  Tp-  16  sec 


Figure  20.  Amplification  factor  sensitivity  to  wave  direction,  Tp=  18  sec 
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Proposed  modifications  were  evaluated  in  terms  of  wave  height  in  the  alterna¬ 
tive  plan  relative  to  wave  height  in  the  existing  harbor  at  three  locations  along  the 
main  wharf  (Basins  1,3,  and  4).  Relative  wave  heights  were  plotted  against  wave 
approach  direction  for  each  harbor  alternative.  The  two  peak  wave  periods 
considered,  10  and  16  sec,  are  representative  of  wind  wave  and  swell  conditions 
at  Kaumalapau  Harbor.  Results  for  the  dogleg  breakwater  extension  are  given  in 
Figures  21  and  22.  The  addition  of  the  dogleg  extension  significantly  reduces 
wave  heights  within  the  harbor.  Waves  with  10-sec  periods  approaching  from  the 
northwest  were  particularly  affected,  resulting  in  wave  heights  43  to  74  percent  of 
those  in  the  existing  harbor.  The  same  wave  condition  approaching  from  the  west 
and  southwest  was  reduced,  but  to  a  lesser  degree,  resulting  in  wave  heights  60  to 
98  percent  those  in  the  existing  harbor.  A  single  case  of  increased  amplification 
occurred  at  the  east  end  of  the  main  wharf  (Basin  4)  for  an  incident  wave  direc¬ 
tion  of  221  deg. 

Longer  period  swell  (Figure  22)  also  showed  a  marked  decrease  in  wave 
height  with  the  addition  of  the  dogleg  extension,  but  with  greater  consistency 
through  all  incident  directions.  For  waves  approaching  between  315  and  201  deg, 
wave  heights  along  the  main  wharf  ranged  between  65  and  80  percent  of  those  in 
the  existing  harbor. 

The  addition  of  wave  absorber  along  the  northeast  or  south  shore  of  the  harbor 
had  little  effect  on  wave  heights  at  the  wharf.  In  all  cases  involving  wave 
absorbers,  the  average  reduction  in  wave  height  along  the  wharf  was  no  greater 
than  10  percent.  Although  wave  absorbers  had  little  impact  on  the  wharf  area, 
they  did  reduce  wave  heights  in  a  very  localized  area  immediately  adjacent  to  the 
absorbing  shore. 

Similar  results  for  the  straight  breakwater  extension  are  shown  in  Figures  23 
and  24.  For  10-sec  waves  approaching  from  the  northwest,  the  addition  of  the 
straight  breakwater  extension  without  further  modification  gives  wave  heights 
along  the  pier  between  45  and  71  percent  of  those  within  the  existing  harbor.  The 
same  waves  approaching  from  the  southwest  show  wave  heights  that  are  between 
61  and  113  percent  those  with  the  existing  harbor,  where  values  greater  than 
100  percent  represent  an  increase  in  wave  amplification.  Longer  period  swell 
shows  less  directional  sensitivity,  giving  a  more  focussed  range  of  wave  heights 
between  69  and  90  percent  of  those  for  existing  harbor  conditions. 

Based  on  results  for  dogleg  modifications,  wave  absorber  was  added  only 
along  the  north  shore  for  straight  breakwater  designs.  Again  the  average  reduc¬ 
tion  in  wave  height  along  the  pier  did  not  exceed  10  percent. 

In  order  to  isolate  effects  of  the  wave  absorber  on  harbor  response,  absorber 
combinations  with  the  existing  breakwater  were  also  evaluated  (Figures  25  and 
26).  Results  show  that  wave  heights  along  the  pier  are  reduced  by  no  more  than 
1 1  percent  for  all  cases  where  wave  absorber  was  present.  Also,  the  addition  of  a 
south  breakwater  had  little  effect  on  wave  heights  along  the  pier,  showing  a 
maximum  reduction  in  height  of  19  percent  with  an  average  reduction  of  approxi¬ 
mately  2  percent. 
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WAVE  HEIGHT  RELATIVE  TO  EXISTING  HARBOR  (%) 


Figure  21.  Wave  height  relative  to  existing  harbor,  dogleg  breakwater 
extension,  Tp  =  1 0  sec 
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WAVE  HEIGHT  RELATIVE  TO  EXISTING  HARBOR  (%) 


Figure  22.  Wave  height  relative  to  existing  harbor,  dogleg  breakwater 
extension,  Tp  =  1 6  sec 
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WAVE  HEIGHT  REU\TIVE  TO  EXISTING  HARBOR  (7.) 
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WAVE  HEIGHT  RELATIVE  TO  EXISTING  HARBOR  (7.) 


Figure  26.  Wave  height  relative  to  existing  harbor,  existing  breakwater,  Tp- 
1 6  sec 
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Amplification  factor  contour  plots  illustrate  the  effect  of  each  of  the  break¬ 
water  alternatives  on  waves  entering  the  harbor  (Figures  27  and  28).  The  effect 
of  the  wave  absorber  along  the  northeast  shore  also  is  illustrated. 

In  summary,  the  numerical  model  simulations  of  harbor  alternatives  indicate 
that  a  45.7-m  (150-ft)  extension  to  the  existing  breakwater  will  reduce  short  wave 
heights  at  the  wharf  to  about  60  to  80  percent  of  those  in  the  existing  harbor.  The 
dogleg  breakwater  extension  is  a  little  more  effective  than  the  straight  extension 
for  reducing  wave  action  at  the  wharf.  Wave  absorbers  along  the  northeast  or 
south  shores  of  the  embayment  have  minimal  impact  on  wave  heights  at  the 
wharf. 
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Figure  27.  Amplification  factor  contours,  251 -deg  incident  direction,  Tp  =  10  sec 
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Figure  28.  Amplification  factor  contours,  251 -deg  incident  direction,  Tp  =  16  sec 
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6  Harbor  Oscillations1 


To  evaluate  harbor  resonance  characteristics,  the  HARBD  numerical  model 
was  run  for  the  existing  harbor  and  the  WES-recommended  plan.  Incident  long¬ 
wave  periods  ranged  from  23  to  513  sec  in  very  fine  increments,  as  discussed  in 
Chapter  4.  These  simulations  were  included  because  modifications  to  a  harbor 
can  potentially  lead  to  increased  operational  problems  due  to  harbor  oscillations. 

Amplification  factors  in  the  existing  harbor  along  the  wharf,  Basins  1  through 
4  (Figure  16),  are  shown  as  a  function  of  wave  frequency  in  Figure  29.  Some 
frequencies  produce  a  strong  resonant  amplification,  with  peak  amplification 
factors  between  approximately  5  and  30.  Many  of  the  same  resonant  frequencies 
appear  at  all  four  basins,  though  the  strength  of  amplification  can  vary  consider¬ 
ably  between  basins.  Similar  resonant  frequencies  are  indicated  at  the  other 
basins  with  few  exceptions  (not  shown). 

Amplification  factors  in  the  WES-recommended  alternative  plan  along  the 
wharf  also  show  strong  resonance  at  some  frequencies  (Figure  30).  The  lowest 
frequency  resonance,  at  0.00615  Hz,  is  similar  to  that  in  the  existing  harbor. 
However,  the  other  resonant  peaks  generally  differ  from  the  existing  harbor 
results  in  both  magnitude  and  frequency.  The  WES-recommended  plan  shows 
fewer,  but  generally  higher,  resonant  peaks  than  the  existing  harbor. 

Amplification  factor  and  phase  contour  plots  for  the  six  highest  resonant  peaks 
show  oscillation  patterns  in  the  existing  harbor.  In  the  amplification  factor  plots, 
areas  of  high  amplification  are  evident  as  darker  shades  of  gray  (Figure  31). 
Corresponding  phase  contours  are  shown  in  Figure  32.  In  a  pure  standing  wave, 
the  phase  values  are  constant  between  nodes  and  they  shift  180  deg  across  a  node. 
Areas  in  which  phase  contours  are  tightly  bunched  indicate  nodal  areas.  Veloci¬ 
ties  and  particle  excursions  at  nodes  are  relatively  large  and  vertical  motions  are 
relatively  small.  Thus,  nodal  lines  in  Figure  32  coincide  with  low  amplification 
factors  in  Figure  31.  The  phase  plots  also  indicate  areas  of  the  harbor  which  rise 
and  fall  together  during  the  resonant  condition  (same  gray  shade).  Thus,  the 
oscillation  patterns  can  be  interpreted. 


’Written  by  Edward  F.  Thompson  and  Lori  L.  Hadley. 
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Figure  31 .  Resonant  long  wave  amplification  factor  contours,  existing  harbor 
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Figure  32.  Resonant  long-wave  phase  contours,  existing  harbor 
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From  the  figures,  it  is  clear  that  peak  ,4  at  161.1  sec  represents  the  Helmholtz 
(or  grave )  mode  of  oscillation,  in  which  the  entire  harbor  rises  and  falls  in  unison. 
The  phase  value  is  constant  over  the  whole  embayment.  The  67. 1-sec  resonant 
period,  peak  B,  represents  a  simple  rocking  oscillation  between  the  north  and 
southeast  areas  of  the  harbor.  A  single  nodal  line  runs  across  the  harbor  in  a 
northeast-southwest  direction.  The  shorter  period  oscillations  are  more  complex 
patterns,  though  they  generally  indicate  a  strong  nodal  area  at  or  near  the  wharf 
and  a  strong  amplification  west  of  the  west  wharf  end  and/or  just  east  of  the  east 
wharf  end.  The  elbow  at  the  intersection  of  the  breakwater  and  the  wharf  and  the 
northeast  corner  of  the  harbor  appear  to  be  natural  antinodes  for  long-wave 
resonance.  Peaks  D  and-  G,  which  are  not  included  in  the  figures,  are  weaker 
resonances,  in  which  the  north  end  of  the  embayment  acts  in  phase  and  there  are 
no  nodal  lines  intersecting  the  wharf.  The  peak  I  resonance  (not  shown)  is  an 
interesting  pattern  similar  to  that  shown  below  for  peak  I  in  the  WES-recom¬ 
mended  plan. 

Amplification  factor  and  phase  contour  plots  for  resonant  peaks  in  the  WES- 
recommended  plan  show  similar  oscillation  patterns  for  peaks  A,  B,  C,  E,  and  F 
(Figures  33  and  34).  The  stronger  amplifications  for  the  peak  B,  C,  and  F  cases 
(relative  to  the  existing  harbor)  can  be  attributed  to  the  greater  constriction  of  the 
harbor  by  the  longer  breakwater  in  the  WES-recommended  plan.  There  is  no 
strong  resonance  in  the  WES-recommended  plan  comparable  to  peak  H  in  the 
existing  harbor,  and  the  pattern  for  peak  /  is  included  instead.  The  peak  I 
resonance  shows  an  interesting  nodal  pattern  which  nearly  intersects  the  western 
part  of  the  wharf. 

The  long-wave  amplification  factors  shown  here  should  be  viewed  as  conserva¬ 
tively  high  for  several  reasons.  Wave  reflection  coefficient  at  all  solid  boundaries 
was  set  to  1.0.  Bottom  friction  was  set  to  0.0.  Also,  the  breakwater  is  repre¬ 
sented  as  a  solid  barrier;  but  for  harbor  oscillation  wave  periods,  significant 
energy  may  be  transmitted  through  it. 

A  more  quantitative  comparison  between  the  existing  harbor  and  the  WES- 
recommended  alternative  plan  can  be  obtained  by  averaging  amplification  factors 
across  the  entire  range  of  long-wave  frequencies.  The  root-mean-square  (RMS) 
amplification  factor  was  computed  for  each  plan  (Table  14).  The  RMS  is  used 
because  squared  amplification  factors  are  indicative  of  wave  energy,  a  more 
relevant  basis  for  comparison  than  wave  height.  The  WES-recommended  plan 
gives  a  higher  RMS  amplification  factor  than  the  existing  harbor  at  all  four  basin 
locations  along  the  wharf.  As  might  be  expected,  the  biggest  increase  between 
plans  (20  percent)  is  at  the  western  basins,  closest  to  the  breakwater.  Differences 
between  the  existing  and  WES-recommended  plans  may  be  exaggerated  because 
of  the  very  high,  narrow  resonant  peaks.  Experience  at  other  harbors  has  shown 
that  such  narrow  peaks  are  often  strongly  attenuated  when  bottom  friction  is  added 
to  the  model.  Repetition  of  the  Kaumalapau  Harbor  runs  with  a  representative 
value  of  bottom  friction  was  beyond  the  scope  of  this  study. 
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Figure  33.  Resonant  long-wave  amplification  factor  contours,  WES-recommended  plan 
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Figure  34.  Resonant  long-wave  phase  contours,  WES-recommended  plan 
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Table  14 

Comparison  of  Long  Wave  Response  in  Existing  and  WES-Recommended 
Harbor  Plans 

Basin 

RMS  Amplification  Factor 

Existing 

Harbor 

WES-Recommended  Plan 

Percent  Difference 
(WES  Plan  -  Existing) 

1  (Mid-Wharf) 

3.13 

3.49 

+12 

2  (W.  Comer) 

4.13 

4.94 

+20 

3  (W.  Wharf) 

3.33 

3.98 

+20 

4  (E.  Wharf) 

3.95 

4.18 

+  6 

Although  the  WES-recommended  plan  is  prone  to  stronger  oscillations  than  the  existing 
harbor,  it  seems  likely  that  the  differences  will  have  little  operational  significance.  Even 
with  the  WES-recommended  plan,  the  harbor  is  quite  open.  Thus  the  level  of  wind  wave 
and  swell  energy  at  the  wharf  will  continue  to  be  fairly  high  (by  typical  harbor  standards). 
Wind  waves  and  swell  can  be  expected  to  continue  as  a  primary  operational  concern  which 
will  probably  outweigh  any  practical  problems  due  to  increased  oscillations. 
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7  Physical  Model  Design1 


Introduction 

An  undistorted,  three-dimensional  physical  model  of  Kaumalapau  Harbor, 
Hawaii,  was  constructed  to  evaluate  harbor  response  to  short-period  waves  for  a 
prebreakwater  configuration,  existing  conditions,  and  two  alternative  breakwater 
configurations.  Only  the  straight  and  dogleg  alternative  configurations  were  built, 
because  the  other  alternative  harbor  plans  studied  in  the  numerical  model  did  not 
perform  satisfactorily.  The  physical  model  then  was  used  to  “fine  tune”  design  the 
two  preferred  plans.  The  physical  model  also  was  used  to  investigate  and  optimize 
stability  of  the  breakwater  configuration  which  provided  the  best  reduction  in  wave 
energy  at  the  pier. 

Model  Design 

Model  scaling 

The  Kaumalapau  physical  model  was  constructed  at  an  undistorted  linear  scale 
of  1:49,  model  to  prototype.  This  scale  was  chosen  based  upon  the  following 
factors: 

a.  Size  of  the  physical  model  which  can  accurately  reproduce  the  necessary 
wave  patterns  within  physical  limitations  imposed  by  basin  size  and 
hardware  constraints. 

b.  Depth  of  water  necessary  to  generate  realistic  waves  without  excessive 
scale  effects  and  the  ability  to  generate  the  size  of  hurricane  waves  needed. 

c.  Available  CORE-LOC™  (hereafter  referred  to  as  core-loc)  and  dolos  sizes. 

d.  Model  construction  costs. 


'written  by  Gordon  S.  Harkins  and  Ernest  R.  Smith. 
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Following  selection  of  the  linear  scale,  the  model  was  designed  and  operated  in 
accordance  with  Froude’s  model  law  (Stevens  et  al.  1942).  The  scale  relations  used 
for  design  and  operation  of  the  model  are  shown  in  Table  15. 


Table  151 

Model-Prototype  Scale  Relations  (1:49  scale) 

Scale  Relations 

Characteristic 

Dimension 

Model:Prototype 

Length 

l 

1  = 

1:49 

Area 

P 

Ctr  ~ 

1:2401 

Volume 

P 

Vr  ~ 

1:117649 

Time 

V* 

t,  = 

1:7 

Velocity 

Vt 

vr  = 

1:7 

Discharge 

P/t 

<7,  = 

1:16807 

1  Dimensions  are  in  terms  of  length  (/)  and  time  (f). 

The  physical  model  covered  approximately  0.3  sq  km  and  extended  vertically  from 
-29  m  (-95  ft)  to  +8  m  (26.2  ft)  mllw  contour  in  the  prototype  (Figure  35).  The 
+8-m  (+26.2-ft)  contour  was  chosen  to  avoid  overtopping  on  to  the  top  of  the 
model.  However,  water  splashed  on  to  the  overbank  under  extreme  hurricane 
conditions. 

Construction 

The  model  was  designed  using  detailed  bathymetric  and  topographic  data  and 
constructed  with  sand  fill  and  a  concrete  cap  (Figure  36).  In  the  vicinity  of  the 
existing  breakwater,  the  model  was  constructed  in  a  pre-breakwater  configuration. 
This  first  configuration  allowed  the  determination  of  the  design  wave  heights  for 
model  armor  unit  selection.  The  existing  and  two  alternative  breakwater  configura¬ 
tions  then  were  built  above  the  pre-breakwater  contours  as  described  below. 

Breakwater  design 

The  proposed  breakwater  rehabilitation  was  designed  to  withstand  storms  from 
the  South  Pacific  with  minimal  repairs.  The  crest  elevation  was  to  be  raised  to 
+6. 1  m  (+20  ft)  mllw  with  a  6. 1-m-wide  (20-ft-wide)  crown  consisting  of  a  cast-in¬ 
place  concrete  ribcap.  The  proposed  length  of  the  structure  was  91.4  m  (300  ft).  It 
was  desired  to  use  existing  breakwater  material  as  a  foundation  and  core  for  the  pro¬ 
posed  breakwater. 


Chapter  7  Physical  Model  Design 


59 


Figure  35.  Physical  model  layout 
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Figure  36.  Photograph  of  physical  model 


Cross  section  design.  The  proposed  breakwater  consisted  of  an  armor  layer 
constructed  of  core-locs,  an  underlayer,  and,  where  necessary,  a  core.  It  was  neces¬ 
sary  in  certain  locations  to  excavate  material  from  the  existing  breakwater  to  place 
the  armor  layer  and  underlayer.  This  excavated  material  was  used  as  a  base  near 
the  end  of  the  proposed  structure.  Core  material  was  added  only  in  locations  whore 
it  was  required  to  raise  the  elevation  of  the  existing  breakwater  to  the  bottom 
elevation  of  the  underlayer. 

The  design  wave  height  specified  by  the  U.S.  Army  Engineer  District  (USAED), 
Honolulu  (1993)  was  a  deepwater  wave  height  H0  of  8.5  m  (28  ft)  or  wave  height  at 
the  structure  Hd  of  7 2  m  (23.6  ft).  Core-loc  size  for  the  armor  layer  was  calculated 
using  the  stability  equation  of  Hudson  (1958): 


W 


Y  aHl 


KAS '-If  cote 


(14) 


in  which  Wa  is  the  weight  of  an  individual  armor  unit,  y0  is  the  specific  weight  of  an 
individual  armor  unit,  Kd  is  the  stability  coefficient,  Sa  is  the  specific  gravity  of  an 
individual  armor  emit  or  stone  relative  to  the  water  in  which  it  is  placed,  Sa  =  y  Jyw 
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in  which  yw  is  the  specific  weight  of  water,  and  dis  the  angle  of  the  structure  slope 
measured  from  horizontal  in  degrees. 

The  water  depth  at  the  breakwater  exceeds  1.5 Hd  (10.8  m,  35.4  ft).  Therefore, 
core-loc  was  placed  according  to  Engineer  Manual  1 1 10-2-2904  (“Design  of 
Breakwaters  and  Jetties”)  (U.S.  Army  Corps  of  Engineers  1986)  to  a  bottom 
elevation  of  -10.8  m  (-35.4  ft)  mllw. 

Straight  alignment.  Figures  37  and  38  show  cross  sections  of  the  proposed 
breakwater  superimposed  with  the  existing  structure  at  selected  stations.  The 
seaside  slope  of  the  proposed  structure  extended  past  the  existing  structure  at 
the  toe  for  stations  seaward  of  Station  0+25.  This  indicates  that  fill  material 
was  required  to  support  the  structure  sea  side  at  these  stations.  However, 
sufficient  material  exists  on  the  structure  lee  side.  To  make  efficient  use  of 
existing  material,  cross  sections  requiring  fill  on  the  sea  side  were  offset  to  the 
lee  side  to  allow  a  minimum  sea-side  toe  shelf  of  1.5  m  (5  ft).  Figure  39  shows 
a  plan  view  of  the  existing  structure  and  points  indicating  the  minimum  center- 
line  location  for  each  station  to  avoid  sea-side  fill  placement.  An  adjusted  center- 
line  was  drawn,  shown  as  a  dashed  line  in  Figure  39,  that  maintained  the  minimum 
toe  shelf  distance  and  provided  a  straight  breakwater  alignment.  Little  existing 
material  is  present  between  Station  2+25  and  the  terminus  of  the  proposed  structure. 
Therefore,  fill,  obtained  from  the  structure  leeside,  was  required  to  provide  a  base 
for  the  breakwater  in  this  region.  A  cross  section  at  Station  2+00  of  the  revised 
straight  breakwater  alignment  is  shown  in  Figure  40,  and  a  plan  view  sketch  is 
shown  in  Figure  41. 


Figure  37.  Cross  section  of  proposed  structure  at  existing  structure  center 
line,  Station  0+00 
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Figure  38.  Cross  section  of  proposed  structure  at  existing  breakwater 
center  line,  Station  2+00 


Figure  39.  Minimum  center  line  offset  distances  and 
revised  centerline  for  straight  breakwater 
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HARBOR  SIDE 


Ribcap 

Vi 


6.1  m 


SEASIDE 


STATION  2  +  00.0 


Proposed  Breakwater 
Existing  Breakwater 
Excavated  Existing  Breakwater 


W,  =  31.4  tonnes 
Wj  =3.1  to  4.8  tonnes 
W3  =  0.5  to  68  kg 
W«  =  4.8  to  10.0  tonnes 


Figure  40.  Cross  section  of  proposed  breakwater  at  revised  center  line, 
Station  2+00 


Dogleg  alignment.  An  alternative  to  the  straight  breakwater  alignment  was  a  dog¬ 
leg  configuration.  The  breakwater  head  position  was  limited  by  a  navigational  buoy 
and  the  harbor  basin.  A  configuration  was  developed  that  provided  the  maximum 
protection  within  the  constraints  of  navigation.  The  dogleg  portion  consisted  of 
1 1.9  m  (39  ft)  of  the  seaward  portion  of  the  structure.  The  dogleg  originated  at 
Station  2+34  and  extended  toward  the  harbor  1 1.9  m  (39  ft)  at  an  angle  of  29  deg 
with  the  straight  alignment  center  line  (45  deg  with  the  original  center  line. 

Figure  42). 


Figure  42.  Sketch  of  dogleg  breakwater 
configuration  (Plan  4) 


Equipment 

Wave  generators 

Two  wavemakers  were  connected  to  provide  a  9.  l-m-(30-ft-)  wide  crest  width. 
This  corresponds  to  the  prototype  width  of  approximately  457.2  m  (1,500  ft),  which 
is  the  same  order  of  magnitude  as  the  harbor  entrance.  The  wavemakers  used  in  this 
study  were  hydraulically  powered,  piston  wavemakers.  They  were  capable  of 
making  waves  from  one  direction  (i.e.,  unidirectional)  at  a  time.  To  generate  waves 
from  different  directions,  the  wavemakers  were  physically  reoriented. 
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Heights  of  waves  generated  from  a  wavemaker  are  a  function  of  the  water  depth 
at  the  wavemaker.  To  generate  waves  of  significant  height  and  to  minimize 
spurious  wave  components,  the  wave  machines  were  located  in  the  prototype  depth 
of  -29  m  (-95  ft). 

The  wavemakers  were  computer  controlled  and  generated  irregular  waves. 
Irregular  waves  have  varying  wave  height  and  period  and  provide  more  realistic  sea 
conditions  than  monochromatic  waves. 

Wave  gauges 

Parallel  wire  capacitance  wave  gauges  were  used  to  measure  and  record  model 
wave  heights.  Output  of  gauges  is  directly  proportional  to  their  submergence  depth 
in  water  and  calibration  curves  for  each  wave  gauge  were  recorded  each  morning, 
by  automatically  stepping  wave  rods  throughout  the  water  column  over  1 1  steps. 
Two  wave  gauge  configurations  were  used  throughout  model  experiments. 

The  first  configuration,  called  the  pre-breakwater  configuration,  is  shown  in 
Figure  43.  This  wave  gauge  configuration  utilized  15  gauges  and  was  used  to 
calibrate  the  waves  and  to  determine  wave  heights  along  the  center  line  of  the  pre¬ 
breakwater  configuration.  Gauges  were  numbered  to  correspond  with  the  numerical 
model  basin  locations  (Figure  16).  Gauge  locations  labeled  1  through  7  correspond 
with  the  basin  locations  labeled  identically.  Gauges  1  through  4  were  the  principal 
gauges  of  interest  and  were  located  along  the  main  pier  with  gauge  2  located  closest 
to  the  breakwater,  and  gauges  3  and  1  in  the  center  of  the  pier,  and  gauge  4  farthest 
from  the  breakwater.  Gauge  5  was  located  in  the  center  of  the  mouth  of  the  bay. 
Gauge  7  was  located  approximately  100  m  (330  ft)  from  the  pier  and  gauge  8  was 
positioned  30  m  (100  ft)  from  the  pier  and  50  m  (160  ft)  from  the  base  of  the  cliff 
wall  to  quantify  reflection  of  wave  energy  off  the  near-vertical  walls.  If  numerical 
model  results  had  indicated  that  wave  absorbers  were  effective,  then  this  gauge 
would  have  been  used  to  determine  changes  between  cases  with  a  wave  absorber 
versus  those  without  a  wave  absorber.  Gauge  9  corresponds  to  basin  location  1 1 
and  initially  was  going  to  be  used  to  measure  incident  wave  height  for  one  of  the 
preliminary  plans  for  a  jetty  located  along  the  south  side  of  Kaumalapau  Harbor. 
Results  from  this  plan  were  not  satisfactory  and  were  not  continued  beyond  numeri¬ 
cal  model  studies.  Gauges  10  through  13  were  located  along  the  center  line  of  the 
existing  breakwater  location.  The  model  was  in  a  pre-breakwater  configuration,  so 
that  the  gauges  were  located  above  an  assumed  bottom  depth  which  would  corre¬ 
spond  to  the  depth  present  before  the  breakwater  was  built.  These  gauges  were 
positioned  over  the  center  line  of  the  breakwater  and  were  used  to  define  maximum 
wave  height  occurring  along  the  breakwater.  Gauges  14  and  15  were  used  to 
calibrate  the  incoming  wave  heights  and  are  located  in  a  depth  of  -29  m 
(-95  ft)  mllw. 

The  second  wave  gauge  layout  was  used  during  studies  of  the  alternative  break¬ 
water  configurations  and  can  be  seen  in  Figure  44.  Gauges  1  through  7  were  not 
moved  from  the  first  wave  gauge  configuration  except  for  the  fifth  gauge.  The  fifth 
gauge  was  located  in  two  adjacent  positions  labeled  5  and  5'.  The  location  of 
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Figure  43.  Wave  gauge  configuration  1  for  studies  of 
pre-breakwater  conditions 


gauge  5  corresponds  to  the  numerical  model  basin  location.  For  ease  in  building 
and  removing  the  straight  and  dogleg  breakwater,  the  fifth  gauge  was  repositioned 
slightly  seaward  of  its  initial  position  and  relabelled  5'  to  distinguish  its  second 
position.  The  fifth  gauge  was  positioned  at  location  5'  whenever  the  straight  or  dog¬ 
leg  extensions  were  in  place.  Gauges  6  through  8  correspond  to  the  numerical 
model  basin  locations.  Gauge  9  was  positioned  seaward  of  the  small  cobble  pocket 
beach.  At  one  point,  developers  were  going  to  build  a  man-made  sand  beach  and 
gauge  9  was  located  in  this  position  to  quantify  the  incident  wave  heights  to  the 
shoreline.  There  was  also  talk  of  building  a  floating  private  marina  in  the  vicinity  of 
gauge  10.  Gauge  1 1  was  located  at  the  numerical  model  basin  location  No.  11. 
Gauge  12  was  located  seaward  of  the  breakwater  to  quantify  incident  wave  heights. 
Gauges  13  and  14  were  used  to  measure  the  incident  wave  heights. 

Computer  support 

Wave  machine  signal  generation,  wave  machine  control,  and  data  collection  were 
performed  on  a  Digital  Equipment  Corporation  MicroVax  computer.  This  process 
consists  of  converting  a  digital  control  signal  to  an  analog  voltage  which  is  sent 
simultaneously  to  the  hydraulic  controllers  of  the  wave  machine.  Wave  gauge 


Chapter  7  Physical  Model  Design 


67 


Figure  44.  Wave  gauge  configuration  2  for  studies  of 
existing  and  alternative  harbor  plans 


calibration  and  preliminary  data  analysis  were  also  conducted  on  the  MicroVax 
computer. 

Wave  absorbing  material 

Rolls  of  fiber  wave  absorber  measuring  1.2  m  (4  ft)  wide  and  0.6  m  (2  ft)  high 
were  positioned  behind  the  wavemakers  to  absorb  waves  generated  from  the  back  of 
the  wave  paddles.  Behind  the  rolls  of  fiber  wave  absorbing  material,  2-kg  (5-lb)  to 
15-kg  (30-lb)  stones  on  a  1:2  (vertical horizontal)  slope  up  to  1.5  m  (5  ft)  in  height 
were  positioned  adjacent  to  the  back  wall.  This  sloping  structure  further  reduced 
wave  energy  that  would  have  reflected  off  the  vertical  basin  walls. 

Fiber  wave  absorber  material  was  positioned  around  the  perimeter  of  the  model 
in  front  of  the  wave  machines  at  the  water  surface  to  absorb  reflections  off  the 
Kaumalapau  model. 
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Guide  vanes 


In  nature,  storms  that  generate  wave  events  are  large  in  size  and  can  propagate 
over  long  distances  affecting  large  portions  of  the  coastline.  Finite  length 
wavemakers,  however,  have  end  effects.  To  minimize  end  effects,  aluminum 
barriers  called  guide  vanes  measuring  3  m  (10  ft)  long  by  1.5  m  (5  ft)  tall  and 
approximately  0.5  cm  (0.25  in.)  wide  were  positioned  at  both  ends  of  the  wave 
machine  to  contain  the  wave  energy  generated  by  the  wavemakers.  These  guide 
vanes  limit  the  wave  diffraction  effects  that  occur  at  the  ends  of  the  wave  paddles 
and  also  contain  the  wave  energy,  providing  a  more  uniform  wave  elevation  incident 
to  the  model  study  area. 
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Still-water  Levels 

Still-water  levels  (swl's)  for  wave  action  models  are  selected  so  that  various 
wave-induced  phenomena  that  are  dependent  on  water  depths  are  accurately 
reproduced  in  the  model.  These  phenomena  include  refraction,  diffraction,  and 
overtopping  of  harbor  structures. 

In  most  cases,  for  the  following  reasons,  it  is  desirable  to  select  a  model  swl 
that  closely  approximates  the  higher  prototype  water  levels: 

a.  The  maximum  amount  of  wave  energy  reaching  a  coastal  area  normally 
occurs  during  the  higher  water  phase  of  the  local  tidal  cycle. 

b.  Most  storms  moving  onshore  are  characteristically  accompanied  by  a 
higher  water  level  due  to  wind,  tide,  and  shoreward  mass  transport. 

c.  The  selection  of  a  high  swl  helps  minimize  model  scale  effects  due  to 
viscous  bottom  friction. 

d.  When  a  high  swl  is  selected,  a  model  investigation  tends  to  yield  more 
conservative  results. 

The  Hawaiian  Islands  experience  two  high  and  two  low  tides  daily.  The  tides 
are  semidiurnal  with  pronounced  diurnal  inequalities.  Tidal  data  representative  at 
the  site  are  shown  in  Table  16  (USAED,  Honolulu  1993): 


Written  by  Gordon  S.  Harkins. 
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Table  16 

Kaumalapau  Tidal  Data 

Tide  Description 

Elevation 

m  (ft) 

Highest  tide  (estimated) 

+1 .06  (+3.5) 

Mean  higher  high  water 

+0.67  (+2.2) 

Mean  high  water 

+0.52  (+1.7) 

Mean  tide  level 

+0.29  (+0.95) 

Mean  low  water 

+0.06  (+0.2) 

Mean  lower  low  water 

0.0  (0.0) 

Lowest  tide  (estimated) 

-0.3  (-1.0) 

Swl's  of  0.67  m  (2.2  ft)  and  1.52  m  (5.0  ft)  mllw  were  selected  by  POD  for 
use  in  experiments  of  Kaumalapau  Harbor.  The  0.67-m  (2.2-ft)  value  was 
representative  of  mean  higher  high  water  (mhhw)  and  was  used  during  operational 
wave  condition  experiments  from  the  various  directions.  The  1.52-m  (5.0-ft) 
level  consisting  of  an  astronomical  tide  of  1.06  m  (3.5  ft),  a  water  level  rise  due 
to  a  reduction  in  atmospheric  pressure  of  0.43  m  (1.4  ft),  and  a  water  level  rise 
due  to  a  storm  surge  of  0.03  m  (0. 1  ft),  was  used  during  experiments  of  severe 
storm  wave  conditions  associated  with  passage  of  the  design  hurricane. 

Factors  Influencing  Selection  of  Experiment  Wave 
Characteristics 

In  planning  the  physical  model  experimental  program  for  Kaumalapau  Harbor, 
it  was  necessary  to  select  heights,  periods,  and  directions  for  the  experimental 
waves  that  would  allow  realistic  study  of  the  two  proposed  alternative  breakwater 
configurations  based  upon  typical  wave  conditions  and  the  design  hurricane  waves 
as  seen  at  the  site.  WIS  data  were  available  to  the  southwest  of  the  site  and 
provided  the  primary  means  of  determining  wave  height  and  period.  Wave 
direction  was  determined  from  the  WIS  data  and  numerical  modeling  results. 

Since  explicit  wave  criteria  for  reduction  in  wave  heights  at  the  wharf  were  not 
known,  typical  waves  found  at  Kaumalapau  were  examined.  Hurricane  waves 
were  used  in  the  breakwater  stability  studies. 

Wave  refraction 

When  waves  move  into  water  of  gradually  decreasing  depth,  transformations 
take  place  in  all  wave  characteristics  except  wave  period  (to  the  first  order  of 
approximation).  The  most  important  transformations  with  respect  to  the  selection 
of  experimental  wave  characteristics  for  the  study  are  the  changes  in  wave  height 
and  direction  of  travel  due  to  the  phenomenon  referred  to  as  wave  refraction. 
When  the  refraction  coefficient  ( Kr )  is  determined,  it  is  multiplied  by  the  shoaling 
coefficient  (Ks),  which  gives  a  conversion  factor  for  transforming  deepwater  wave 
heights  to  shallow-water  values.  The  shoaling  coefficient,  a  function  of  wave 
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length  and  water  depth,  can  be  obtained  from  the  Shore  Protection  Manual 
(1984).  The  change  in  wave  height  and  direction  may  be  determined  by  using  the 
numerical  Regional  Coastal  Processes  Wave  Transformation  Model  (RCPWAVE) 
developed  by  Ebersole,  Cialone,  and  Prater  (1986). 

Due  to  the  very  deep  water  levels  just  offshore  of  Kaumalapau  Harbor,  a  wave 
refraction  analysis  was  not  deemed  essential.  Deepwater  waves  generally  begin 
refracting  at  a  depth  of  about  one  half  their  wave  length  (deepwater  wave  length 
La  is  defined  as  gf!2  ).  A  wide  range  of  wave  conditions  was  studied  in  the 
model.  Waves  were  generated  in  the  -29-m  (-95-ft)  model  pit.  From  this  point, 
the  model  contours  refracted  the  wave  trains  to  the  shore. 

Wave  and  storm  data  selected 

Prevailing  winds  in  the  Hawaiian  Islands  are  the  northeasterly  trade  winds, 
which  occur  approximately  90  to  95  percent  of  the  time  during  the  summer 
months  (May-October),  and  55  to  65  percent  of  the  time  between  November  and 
April,  with  speeds  of  16  to  32  km/hr  (10  to  20  mph).  Storm  conditions  generally 
result  when  a  breakdown  of  the  trade  wind  circulation  occurs,  which  is  relatively 
infrequent. 

Three  classes  of  disturbances  produce  major  storms  in  Hawaii:  cold  fronts, 
low-pressure  passages,  and  true  tropical  storms  or  hurricanes.  Cold  fronts,  which 
occur  during  the  winter  months,  cause  spotty  rainfall  and  gusty  winds.  Low- 
pressure  passages  result  in  heavy  rain,  sometimes  with  strong  winds.  A  low- 
pressure  storm  type  known  as  a  "kona"  storm  usually  occurs  during  the  winter 
months,  and  is  associated  with  strong  and  persistent  southerly  winds  and  intense 
rainfall  on  the  south  and  western  shores  of  the  islands.  Hurricanes,  classified  as 
storms  with  wind  speeds  greater  than  119  km/hr  (74  mph),  are  infrequent,  but 
historical  records  reveal  that  several  have  passed  over  the  island  of  Lanai. 

Measured  prototype  wave  data  covering  a  sufficiently  long  duration  from  which 
to  base  a  comprehensive  statistical  analysis  of  deepwater  wave  conditions  for  the 
Kaumalapau  Harbor  area  were  not  available.  However,  statistical  wave  hindcast 
estimates  representative  of  this  area  were  obtained  from  the  WIS  studies,  which 
include  a  20-year  hindcast  period  (1956  to  1975).  The  hindcast  (Corson  et  al. 
1986)  was  obtained  at  WIS  Station  34  (20.00  deg  north,  160.00  deg  west)  and 
Station  35  (17.99  deg  north,  157.90  deg  west)  in  the  North  Pacific  Ocean  south¬ 
west  of  the  island  of  Lanai.  The  WIS  data  do  not  include  waves  generated  by 
hurricanes  or  southern  hemisphere  swell  conditions. 

Hurricane  wave  conditions  were  provided  in  Kaumalapau  Harbor  Navigation 
Improvements  Island  of  Lanai,  Hawaii  -  Reconnaissance  Report  (USAED, 
Honolulu  1993),  which  defined  the  design  conditions  as  a  deepwater  significant 
wave  height  of  8.5  m  (27.9  ft)  with  a  9.8-sec  period  based  upon  numerical  model 
predictions.  This  was  assuming  a  worst  case  Hawaiian  hurricane  approaching 
from  the  south  to  southwest  direction. 
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Wave  selection 


The  unidirectional  wave  machines  were  positioned  in  three  directions  to 
reproduce  and  encompass  the  typical  wave  climate  at  Kaumalapau  Harbor. 
Results  from  the  WIS  data  and  model  results  provided  by  HARBD  on  the 
sensitivity  of  the  harbor  to  wave  direction  were  used  to  determine  the  wave 
maker  positioning.  Waves  approaching  from  221,  251,  and  291  deg  referenced 
to  north  were  chosen.  Figure  45  shows  the  wavemaker  positions  and  the  wave 
directions  chosen.  The  wave  period  and  wave  height  for  each  direction  were 
chosen  through  a  review  of  the  data  from  the  WIS  Stations  34  and  35  described 
earlier.  The  combined  number  of  occurrences  from  the  two  stations  bracketed 
by  direction  can  be  seen  in  Table  17.  The  total  number  of  occurrences  from  both 
WIS  stations  is  116,880. 


Figure  45.  Incident  wave  directions  for  physical  model  studies 
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Using  the  WIS  data  above,  POD  agreed  to  the  wave  conditions  shown  in 
Table  18  for  the  harbor  wave  action  portion  of  the  study.  Incident  direction  and 
significant  wave  height  (H^  )are  defined  at  -29  m  (-95  ft)  mllw.  These  waves 
were  reproduced  at  the  0.67-m  (2.2-ft)  swl. 


Table  18 

Physical  Model  Wave  Conditions  Chosen 

Wave  Direction1 
deg 

Wave  Period 

sec 

Hm0Wave  Height1 
m  (ft) 

221 

8.0 

2.5  (8.2), 3.5  (11.5),  4.5  (14.8) 

10.0 

2.5  (8.2), 5.5  (18.0) 

12.0 

1 .5  (4.9),  2.5  (8.2),  3.5  (1 1 .5) 

14.0 

1.5  (4.9),  2.5  (8.2),  3.5  (11.5) 

16.0 

1.5  (4.9),  2.5  (8.2),  3.5  (11.5) 

20.0 

1.5  (4.9),  2.5  (8.2) 

251 

8.0 

2.5  (8.2),  3.5  (11.5),  4.5  (14.8) 

10.0 

2.5  (8.2),  4.5  (14.8),  5.0  (16.4),  5.5 
(18.0),  6.0  (19.7) 

12.0 

2.5  (8.2),  3.5  (11.5) 

14.0 

3.5  (11.5) 

16.0 

3.5  (11.5) 

20.0 

1 .5  (4.9),  2.5  (8.2) 

291 

8.0 

2.5  (8.2),  3.5  (11.5),  4.5  (14.8) 

10.0 

2.5  (8.2),  3.5  (11.5) 

12.0 

2.5  (8.2),  3.5(11.5),  4.5  (14.8) 

14.0 

2.5  (8.2),  3.5  (1 1 .5),  4.5  (1 4.8),  5.5 
(18.0) _ 

1  Defined  at  the  -29-m  (-95-ft)  contour 

Ocean  waves  consist  of  waves  of  more  than  one  height  and  period.  To  reproduce 
typical  ocean  waves,  spectral  seas  were  reproduced  in  the  physical  model.  A 
spectral  parameter  called  the  groupiness  parameter  was  adjusted  to  simulate  how  far 
the  waves  had  traveled  from  the  generating  source.  The  groupiness  of  waves  is 
important  for  long-period  swell  waves  (wave  periods  from  12  to  25  sec)  that  have 
traveled  from  distant  generating  sources.  This  type  of  phenomenon  is  important  for 
both  southern  and  northern  swell. 

The  signals  used  to  control  the  wave  paddles  were  calibrated  by  comparing  the 
target  Hm0  wave  height  versus  the  measured  Hm0  wave  height  at  gauges  14  and  15 
(Figure  43).  A  gain  factor  was  applied  to  the  spectra  and  the  waves  were  rerun  and 
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reanalyzed.  Table  19  shows  the  target  and  reproduced  wave  height  comparisons  for 
a  selection  of  experimental  wave  conditions. 


Table  19 

Target  Vesus  Measured  Wave  Heights  for  Waves  from 

251  deg 

Wave  Period 

sec 

Hm0  Wave  Height 

mm _ 

Target 

Measured 

8 

2.5  (8.2) 

2.5  (8.2) 

8 

3.5  (11.5) 

3.4  (11.2) 

8 

4.5  (14.8) 

4.4  (14.4) 

10 

2.5  (8.2) 

2.5  (8.2) 

10 

6.0(19.7) 

6.0  (19.7) 

10 

5.5  (18.0) 

5.5  (18.0) 

10 

5.0(16.4) 

4.9(16.1) 

10 

4.5  (14.8) 

4.4  (14.4) 

12 

2.5  (8.2) 

2.5  (8.2) 

12 

3.5  (11.5) 

3.5  (1 1 .5) 

14 

3.5  (11.5) 

3.3  (10.8) 

16 

3.5  (1 1 .5) 

3.2  (10.5) 

20 

1.5  (4.9) 

1.5  (4.9) 

20 

2  5  18  71 

2.5  18  71 

Harbor  Layout  Configurations  Studied 

Four  harbor  layout  configurations  were  studied  in  the  physical  model.  Plan  1  was 
the  pre-breakwater  condition,  which  was  used  to  determine  appropriate  core-loc 
armor  unit  sizes  based  upon  the  measured  wave  height  along  the  assumed  existing 
breakwater  center  line.  Wave  response  at  the  pier  area  was  analyzed  for  the 
remaining  three  breakwater  configurations.  Plans  2-4  shown  in  Figure  46.  The 
figure  shows  the  +3-m  (10-ft)  contour  of  the  existing  structure,  the  +6-m  (20-ft) 
contour  of  the  straight  and  dogleg  breakwater  and  the  -9-m  (30-ft)  contour  of  each 
structure.  The  existing  breakwater  condition,  Plan  2,  was  evaluated  first  to  provide 
a  baseline  against  which  the  alternative  configurations  could  be  compared.  The 
second  breakwater  configuration  studied  was  called  the  straight  breakwater  (Plan  3). 
This  was  referred  to  as  the  pre-Ewa  breakwater  with  a  length  of  122  m  (400  ft)  as 
shown  in  Figure  41.  The  last  alternative  studied  was  termed  the  dogleg  breakwater. 
Plan  4.  For  Plan  4,  the  last  15  m  (49.2  ft)  of  the  breakwater  center  line  was  aligned 
30  deg  towards  the  inside  of  the  harbor  as  measured  from  the  center  line  of  the 
straight  breakwater  as  shown  in  Figure  42. 
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Figure  46.  Plan  view  of  breakwater  crest  and  -9.1  m  depth  contour  of  existing,  straight  and 
dogleg  breakwater 
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and  Results1 


Experiment  Series 


Based  upon  the  numerical  model  results,  three  wave  directions  were 
simulated  in  the  physical  model.  The  first  wave  direction  analyzed  was  for  waves 
approaching  the  harbor  directly,  25 1  deg.  This  direction  was  thought  to  provide 
the  highest  waves  that  would  impact  the  breakwater  and  was  used  for  the  pre¬ 
breakwater  experiment  series.  Results  from  this  direction  were  thought  to  provide  a 
good  indication  of  the  merits  of  the  two  proposed  breakwater  configurations.  The 
wave  machines  were  moved  to  the  southerly  direction  and  experiments  on  existing 
conditions  (Plan  2),  straight  (Plan  3),  and  dogleg  (Plan  4)  breakwaters  were 
performed.  Experiments  were  conducted  only  for  existing  conditions  and  the  dog¬ 
leg  breakwater  plan  for  waves  from  291  deg.  Results  from  the  other  two  directions 
indicated  that  the  dogleg  breakwater  was  a  better  alternative  than  the  straight 
breakwater. 

Waves  from  251  deg 

Breakwater  results 

Results  for  the  existing,  straight,  and  dogleg  breakwater  configurations  are  given 
in  Tables  20, 21,  and  22,  respectively.  The  results  show  a  minimal  reduction  in 
deepwater  wave  height,  especially  for  the  longer  period  wind  wave  conditions. 

Results  for  the  8-,  12-,  and  20-sec,  2.5-m  (8.2-ft)  waves  are  shown  graphically 
in  Figure  47.  The  gauges  are  arranged  in  order  of  increasing  distance  from  the 
breakwater.  Results  for  the  straight  or  dogleg  breakwater  do  not  satisfy  the  design 
criteria  established  by  POD  for  all  wave  conditions.  The  criteria  was  a  4.7-m 
(15-fit)  offshore  wave  height  reduced  to  1.5  m  (4.9  ft)  inside  the  harbor.  The 
existing  breakwater  condition  reduces  a  2.5-m  (8.2-ft)  offshore  wave  height  to  less 
than  1.5  m  (4.9  ft)  for  only  the  8-sec  wave  period.  For  the  12-  and  20-sec  waves. 
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Table  20 

Results  for  Waves  from  251  deq.  Existing  Breakwater 

Hmu  Wave  Heights,  m  (ft),  at  Indicated  Gauge  Location 

Gauge 

14 

2.334 

(7.657) 

3.24 

(10.63) 

415 

(13.61) 

2.45 

(8.028) 

4.38 

(14.37) 

5.23 

(17.17) 

5.42 

(17.79) 

5.96 

(19.56) 

2.540 

(8.332) 

3.52 

(11.54) 

3.36 
(1 1 .02) 

3.22 

(10.57) 

1.548 

(5.077) 

2.563 

(8.409) 

Gauge 

13 

2.615 

(8.578) 

3.64 

(11.93) 

4.71 

(15.46) 

2.498 

(8.193) 

4.42 

(14.51) 

4.96 

(16.28) 

5.48 

(17.98) 

5.98 

(19.63) 

2.530 

(8.300) 

3.58 

(11.74) 

3.25 

(10.67) 

3.13 

(10.26) 

1.459 

(4.786) 

2.417 

(7.928) 

Gauge 

12 

2.624 

(8.606) 

3.52 

(11.55) 

4.44 

(14.58) 

2.376 

(7.795) 

4.09 

(13.42) 

4.53 

(14.87) 

O)  CD 

5.42 

(17.78) 

2.227 

(7.304) 

3.16 

(10.36) 

3.15 

(10.34) 

2.976 

(9.762) 

1.229 

(4.033) 

2.045 

(6.709) 

Gauge 

11 

2.619 

(8.590) 

3.58 

(11.74) 

4.46 

(14.62) 

2.549 

(8.361) 

4.57 

(15.00) 

4.99 

(16.37) 

5.38 

(17.65) 

5.86 

(19.22) 

2.309 

(7.574) 

3.18 

(10.44) 

2.960 

(9.710) 

2.871 

(9.419) 

1.184 

(3.883) 

1.926 

(6.318) 

Gauge 

10 

1.559 

(5.113) 

2.245 

(7.365) 

2.837 

(9.306) 

1.457 

(4.778) 

2.930 

(9.611) 

3.31 

(10.87) 

3.72 

(12.19) 

3.90 

(12.81) 

1.449 

(4.753) 

1.997 

(6.551) 

2.204 

(7.230) 

2.263 

(7.424) 

1.018 

(3.341) 

1.635 

(5.363) 

Gauge 

9 

1.956 

(6.416) 

2.837 

(9.307) 

3.37 

(11.07) 

2.584 

(8.476) 

4.19 

(13.74) 

4.33 

(14.20) 

4.48 

(14.70) 

4.76 

(15.61) 

2.204 

(7.228) 

3.14 

(10.30) 

3.60 

(11.81) 

4.38 

(14.38) 

1.911 

(6.270) 

3.13 

(10.27) 

Gauge 

8 

1.733 

(5.683) 

2.435 

(7.986) 

2.974 

(9.754) 

2.060 

(6.757) 

CD 
CO  ^ 
to  ■»- 
co  C* 

3.93 

(12.89) 

4.21 

(13.80) 

4.52 

(14.84) 

2.57 

(8.421) 

3.54 
(1 1 .62) 

3.20 

(10.50) 

3.92 

(12.85) 

1.380 

(4.527) 

2.276 

(7.465) 

Gauge 

7 

1.510 

(4.953) 

2.076 

(6.811) 

2.562 

(8.404) 

1.705 

(5.593) 

2.784 

(9.131) 

3.03 

(9.95) 

3.23 

(10.61) 

3.28 

(10.76) 

1.973 

(6.472) 

2.541 

(8.334) 

2.245 

(7.366) 

2.410 

(7.904) 

1.157 

(3.794) 

1.905 

(6.249) 

Gauge 

6 

2.513 

(8.242) 

3.45 

(11.30) 

4.33 

(14.22) 

2.768 

(9.079) 

4.97 

(16.30) 

5.54 

(18.16) 

5.98 

(19.61) 

6.65 
(21 .83) 

2.821 

(9.254) 

3.838 

(12.59) 

3.55 

(11.65) 

3.32 

(10.89) 

1.733 

(5.686) 

2.868 

(9.407) 

Gauge 

S 

2.279 

(7.477) 

3.29 

(10.80) 

4.26 

(13.96) 

2.504 

(8.213) 

4.66 

(15.29) 

5.35 

(17.56) 

5.91 

(19.40) 

6.45 

(21.17) 

2.640 

(8.658) 

3.68 

(12.09) 

3.50 

(11.49) 

3.73 

(12.25) 

1.492 

(4.894) 

2.500 

(8.199) 

Gauge 

4 

1.385 

(4.543) 

1.790 

(5.871) 

2.136 

(7.008) 

2.162 

(7.093) 

3.09 

(10.13) 

3.27 

(10.72) 

3.44 

(11.30) 

3.46 

(11.36) 

2.447 

(8.025) 

3.043 

(9.983) 

3.50 

(11.49) 

3.53 

(11.57) 

1.356 

(4.448) 

2.000 

(6.563) 

Gauge 

3 

1.295 

(4.248) 

1.618 

(5.307) 

1.918 

(6.290) 

1.787 

(5.863) 

2.650 

(8.693) 

2.90 

(9.51) 

3.12 

(10.22) 

3.32 

(10.90) 

1.896 

(6.220) 

2.475 

(8.117) 

3.25 

(10.65) 

2.775 

(9.102) 

1.579 

(5.179) 

2.327 

(7.632) 

Gauge 

2 

0.686 

(2.249) 

0.9 

(2.952) 

1.108 

(3.635) 

0.891 

(2.922) 

1.389 

(4.557) 

1.54 

(5.04) 

1.69 

(5.54) 

1.84 

(6.05) 

1.188 

(3.897) 

1.539 

(5.048) 

1.697 

(5.566) 

2.015 

(6.611) 

1.262 

(4.142) 

1.882 

(6.173) 

Gauge 

1 

1.447 

(4.746) 

1.881 

(6.171) 

2.213 

(7.259) 

2.129 

(6.982) 

3.113 

(10.21) 

3.320 

(10.89) 

3.543 

(11.62) 

3.640 

(11.94) 

2.481 

(8.138) 

3.128 

(10.26) 

2.904 

(9.525) 

3.210 

(10.53) 

2.180 

(7.151) 

3.317 

(10.88) 

C 

O 

♦3 

=5 

c 

0 

Height 
m  (ft) 

2.50 

(8.2) 

3.50 

(11.48) 

4.50 

(14.76) 

2.50 

(8.2) 

4.39 

(14.4) 

4.94 

(16.20) 

5.46 

(17.90) 

5.97 

(19.60) 

2.50 

(8.2) 

3.50. 

(1 1 .48) 

3.50 
(1 1 .48) 

3.50 
(1 1 .48) 

S  5 

2.50 

(8.2) 

o 

0) 

> 

re 

2 

Period 

sec 

CO 

o 

CM 

•'t 

CD 

20 

Chapter  9  Physical  Model  Experiments  and  Results 


79 


Table  21 

Results  for  Waves  from  251  deg,  Straight  Breakwater 


Gauge 

14 

2.36 

(7.73) 

3.26 

(10.69) 

4.14 

(13.57) 

2.46 

(8.08) 

4.38 

(14.36) 

4.91 

(16.09) 

5.40 

(17.70) 

5.90 

(19.36) 

2.52 

(8.27) 

3.48 

(11.40) 

3.27 

(10.73) 

3.12 

(10.24) 

1.56 

(5.13) 

2.62 

(8.58) 

Gauge 

13 

2.76 

(9.02) 

3.80 
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4.81 

(15.93) 

2.57 

(8.44) 

4.45 

(14.59) 

5.02 

(16.46) 

5.48 

(17.99) 

5.90 

(19.59) 

2.50 

(8.19) 

3.51 
(1 1 .50) 

3.23 
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1.41 
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2.39 

(7.85) 
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5.33 

(17.49) 
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3.31 

(10.87) 

3.34 

(10.97) 

2.95 

(9.66) 

1.08 

(3.53) 

1.78 

(5.85) 

Gauge 
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Table  22 

Results  for  Waves  from  251  deg.  Dogleg  Breakwater  _ 

Wave  Condition  Hm0  Wave  Heights,  m  (ft),  at  Indicated  Gauge  Location 
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Comparison  of  Harbor  Layouts  for  Waves  from  251  deg 

Peak  Period  =  8  sec,  HmO  =  2.5  m  _ 
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3  1 

Gage  Number 


Comparison  of  Harbor  Layouts  for  Waves  from  251  deg 

Peak  Period  =  12  sec,  HmO  =  2.5  m 


-m-  Existing  Condition 
-o-  305  m  Straight  Extension 
-a-  305  m  Dog-leg  Extension 


Figure  47.  Comparison  of  existing  conditions  with  the  two  alternative  layout  configurations  - 
waves  from  251  deg 
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the  wave  height  is  only  less  than  1.5  m  (4.9  ft)  at  gauge  2  for  the  12-sec  wave 
period.  Longer-period  waves  are  only  minimally  affected  by  the  straight  or  dogleg 
breakwater.  However,  results  are  significantly  better  than  for  existing  conditions. 
One  might  notice  the  high  wave  height  at  gauge  1  for  the  20-sec  wave  condition, 
which  is  higher  than  the  deepwater  incident  condition.  This  appears  to  be  due  to  the 
close  proximity  of  the  gauge  to  the  near  vertical  wall  below  the  pier  area.  Near  a 
purely  reflective  structure,  the  wave  height  can  be  twice  as  high  as  the  incident  wave 
height.  For  shorter  period  waves,  this  phenomenon  was  not  present,  indicating  that 
these  waves  were  not  affected  by  the  wall  and  the  distance  to  the  wall  relative  to  the 
wavelength  was  long  enough  that  doubling  of  the  wave  height  was  not  realized.  It 
also  might  indicate  that  a  standing  wave  pattern  is  being  created.  For  the  higher 
wave  heights,  the  results  at  the  pier  are  correspondingly  higher. 

The  differences  between  the  straight  and  dogleg  extensions  become  more 
obvious  for  the  longer  period  waves.  For  the  20-sec  waves,  the  wave  height  at  the 
pier  for  the  dogleg  breakwater  is  on  average  12  percent  less  than  the  straight 
breakwater  with  an  18-percent  maximum  difference  and  a  5-percent  minimum 
difference  for  gauges  1  through  4. 

Waves  from  221  deg 

The  wave  machines  were  reoriented  to  the  22 1 -deg  direction.  Results  for  the 
existing,  straight,  and  dogleg  configurations  are  shown  in  Tables  23, 24,  and  25, 
respectively.  Results  for  8-,  12-,  and  20-sec  waves  are  shown  in  Figure  48.  The 
shorter,  8-sec  waves  with  a  2.5-m  (8.2-ft)  deepwater  wave  height  are  reduced  on  the 
average  by  a  factor  of  two.  The  deepwater  wave  height  is  reduced  by  a  factor  of  five 
at  gauge  2.  Although  gauge  4  is  located  farthest  from  the  breakwater,  the  maximum 
wave  height  occurs  at  gauge  1  for  12-  and  20-sec  waves.  Refraction  of  the  waves  as 
they  enter  the  harbor  and  reflection  from  the  sidewalls  under  the  pier  is  greater  at 
gauge  1.  For  the  8-  and  12-sec  waves,  both  proposed  plans  show  improvements 
over  existing  conditions.  However,  results  with  20-sec  waves  differ  little  between 
the  three  plans.  The  design  criteria  of  a  1 .5-m  ( 4.9-ft)  wave  or  lower  at  the  pier  is 
not  realized  for  waves  from  221  deg.  All  plans  provide  only  limited  sheltering  for 
waves  from  the  221 -deg  wave  direction  for  the  12-  and  20-sec  waves. 

Waves  from  291  deg 

The  wave  machines  were  moved  for  simulation  of  waves  from  291  deg.  Results 
for  the  existing  conditions  and  the  dogleg  configuration  are  given  in  Tables  26  and 
27,  respectively.  For  waves  from  29 1  deg,  the  straight  breakwater  had  been 
eliminated  as  a  possible  alternative  and  experiments  on  this  configuration  were  not 
conducted.  Unlike  the  first  two  directions,  the  20-sec  wave  was  not  a  typical  wave 
condition  from  this  direction  and  thus  a  shorter,  14-sec  wave  was  generated,  and 
results  are  shown  in  Figure  49.  The  pier  area  is  sheltered  by  both  the  existing  and 
dogleg  breakwater  for  waves  from  291  deg.  An  increased  reduction  in  wave  height 
of  the  dogleg  breakwater  over  existing  conditions  at  gauges  3  and  4  is  shown  in 
Tables  26  and  27.  However,  there  is  not  a  substantial  difference  at  gauge  2. 
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Table  23 

Results  for  Waves  from  221  deg,  Existing  Breakwater 
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Table  24 

Results  for  Waves  from  221  deq,  Straiqht  Breakwater 

Hm0  Wave  Heights,  m  (ft),  at  Indicated  Gauge  Location 

Gauge 

14 

2.24 

(7.24) 

3.12 

(10.23) 

3.98 

(13.05) 

2.12 

(6.96) 

4.81 

(15.79) 

1.34 

(4.39) 

2.23 

(7.32) 

3.14 
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Table  25  "  - - 

Results  for  Waves  from  221  deg.  Dogleg  Breakwater 

Wave  Condition _ _  '  _  Hm0  Wave  Heights,  m  (ft),  at  Indicated  Gauge  Location 

Period  Height  Gauge  Gauge  Gauge  Gauge  Gauge  Gauge  Gauge  Gauge  Gauge  Gaut 
sec  m  (ft)  1  2  3  4  5  6  7  8  9  lo' 
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(4.92)  (5.00)  (2.42)  (3.95)  (3.82)  (4.36)  (4.40)  (3.17)  (3.82)  (4.25)  (2.16)  (3.63)  (3.53)  (4.09)  (4.65) 

2.50  2.27  1.11  1.80  1.90  2.24  2.18  1.52  1.83  2.14  1.09  1.78  1.68  1.98  2.36 

(8.20)  (7.43)  (3.64)  (5.91)  (6.22)  (7.35)  (7.15)  (4.99)  (6.01)  (7.02)  (3.59)  (5.85)  (5.50)  (6.48)  (7.73) 


Table  26 

Results  for  Waves  from  291  deg.  Existing  Breakwater  _ 

Wave  Condition _ _ _ _  Hm0  Wave  Heights,  m  (ft),  at  Indicated  Gauge  Location 
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Table  27  - - — 

Results  for  Waves  from  291  deg,  Dogleg  Breakwater  _ 

Wave  Condition _ _ _ _ _ _ Hm0  Wave  Heights,  m  (ft),  at  Indicated  Gauge  Location 
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Comparison  of  Harbor  Layouts  for  Waves  from  291  deg 

Peak  Period  =  8  sec,  HmO  =  2.5  m 


Gage  Number 


Comparison  of  Harbor  Layouts  for  Waves  from  291  deg 


Peak  Period  =  14  sec,  HmO  =  2.5  m 


Figure  49.  Comparison  of  existing  conditions  with  the  two  alternative  layout 
configurations  for  waves  from  291  deg 
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Summary  of  Results  from  the  Three  Wave 
Directions 


Figures  47  through  49  show  the  importance  of  wave  direction,  period,  distance 
from  the  breakwater,  and  breakwater  configuration  on  measured  wave  heights  along 
the  pier  for  a  2.5-m  (8.2-fit)  deepwater  incident  wave.  Wave  heights  at  the  pier  are 
less  than  1.5  m  (4.9  ft)  for  an  incident  wave  period  of  8  sec  and  a  deepwater  wave 
height  of  2.5  m  (8.2  ft)  for  the  three  directions  studied  and  the  three  breakwater 
configurations  constructed.  From  29 1  deg,  the  wave  height  along  the  pier  also  was 
less  than  1.5  m  (4.9  ft)  for  wave  periods  of  12  and  14  sec  for  the  existing  and 
dogleg  breakwaters  for  an  incident  wave  height  of  2.5  m  (8.2  ft).  The  dogleg  break¬ 
water  provided  increased  reduction  of  wave  heights.  For  a  2.5-m  (8.2-fit)  wave  from 
25 1  deg,  the  wave  height  along  the  pier  was  less  than  1 .5  m  (4.9  ft)  for  12-sec 
waves  for  the  straight  and  dogleg  breakwaters,  but  not  for  existing  conditions.  For 
12-sec  waves  from  221  deg,  Hm0  was  less  than  1.5  m  (4.9  ft)  only  at  gauges  2  and  3. 
Only  gauge  2  showed  wave  heights  less  than  1.5  m  (4.9  ft)  for  20-sec  waves  from 
25 1  and  22 1  deg). 

Obviously,  the  criteria  established  by  POD  for  a  4.7-m  (15.5-ft)  deepwater  wave 
height  not  to  exceed  1.5  m  (4.9  ft)  at  the  pier  were  not  met  by  the  alternatives 
studied.  However,  the  results  indicate  that  both  the  straight  and  dogleg  extensions 
decrease  the  wave  height  at  the  pier  compared  to  existing  conditions.  The  dogleg 
extension  provides  slightly  better  results  over  the  straight  breakwater,  especially  for 
the  longer  period  waves  from  251  deg. 


Effectiveness  of  the  Dogleg  Breakwater  Over  Exist¬ 
ing  Conditions  and  the  Straight  Breakwater  Alter¬ 
native 


To  summarize  the  effectiveness  of  the  dogleg  breakwater,  two  normalizing 
techniques  were  used.  The  first  and  simplest  technique,  termed  the  transmission 
coefficient,  normalizes  the  wave  heights  at  the  pier  by  the  deepwater  wave  height  or 
simply 


(15) 


where  Ht  is  the  wave  height  at  the  i*  gauge  location  and  H0  is  the  average  wave 
height  at  gauges  13  and  14,  the  deepwater  wave  height.  This  technique  is  used 
when  showing  the  results  for  a  single  plan.  If  the  result  is  less  than  unity,  then  there 
is  a  reduction  in  wave  height  at  the  pier  in  comparison  with  the  deepwater  wave 
height.  If  the  wave  height  is  amplified  at  the  pier,  the  value  would  be  greater  than 
unity.  To  simplify  the  results,  all  runs  for  a  particular  period  were  averaged  and  one 
result  for  each  wave  period  is  shown.  The  results  for  the  dogleg  breakwater  are 
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shown  in  Figures  50  through  52  for  the  three  incident  wave  directions.  The  graphs 
show  the  transmission  coefficient  of  the  deepwater  wave  height  for  gauges  located 
along  the  pier  area  (see  Figure  44). 


Figure  50.  Period-averaged  deepwater  wave  height  transmission  coefficient  for 
waves  from  251  deg  with  dogleg  breakwater 


Figure  51 .  Period-averaged  deepwater  wave  height  transmission  coefficient  for 
waves  from  221  deg  with  dogleg  breakwater 
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Figure  52.  Period-averaged  deepwater  wave  height  transmission  coefficient  for 
waves  from  291  deg  with  dogleg  breakwater 


Model  results,  showing  the  transmission  coefficient,  are  summarized  in 
Table  28.  The  table  is  divided  into  the  three  directions  with  typical  wave  periods 
associated  with  each  direction.  The  transmission  coefficient  and  maximum  wave 
height  at  the  29-m  (95-ft)  depth  that  produces  a  1.5-m  (4.9-ft)  average  wave  height 
along  the  pier  is  also  given.  Averaging  was  used  to  summarize  the  results  for 
gauges  1-4  for  varying  wave  heights  for  a  particular  wave  period  and  direction.  For 
example,  for  existing  conditions  for  an  8-sec  wave  period  for  waves  from  221  deg 
the  transmission  coefficients  are  shown  in  Table  29. 

The  maximum  wave  height.  Table  28,  at  the  29-m  (95-ft)  depth  was  calculated 
by  dividing  1.5  m  (4.9  ft)  by  the  transmission  coefficient.  This  number  provides  the 
shipping  companies  an  idea  of  the  maximum  deepwater  wave  height  for  a  given 
wave  direction  and  period  that  will  provide  an  average  1.5-m  (4.9  -ft)  wave  height 
along  the  pier.  Another  useful  way  to  use  this  chart  is  to  calculate  the  wave  height 
along  the  pier  given  incident  wave  height,  period,  and  direction  from  which  waves 
are  approaching  the  harbor  in  deep  water.  For  existing  conditions,  given  a  6-m 
(19.7-ft)  wave  height,  8-sec  period  from  221  deg,  the  average  wave  height  along  the 
pier  is  3  m  (9.8  ft),  conditions  that  would  be  too  rough  for  offloading.  The  trans¬ 
mission  coefficient  for  the  dogleg  configuration  is  always  less  than  the  transmission 
coefficient  produced  by  existing  conditions;  thus,  the  dogleg  design  will  allow  the 
harbor  to  be  used  during  more  energetic  wave  events. 

One  might  question  the  validity  of  averaging  the  results  over  different  simulated 
wave  heights.  Transmission  coefficients  for  waves  with  a  wave  period  of  8  sec  and 
from  a  principal  wave  direction  of  25 1  deg  are  shown  in  Figure  53  for  deepwater 
wave  heights  of  2.5  (8.2  ft),  3.5  (1 1.5  ft)  and  4.5  m  (14.8  ft).  Although  the  results 
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Table  28 

Average  Transmission  Coefficient  and  Maximum  Deepwater  Wave 
Height  That  Provide  for  a  1.5-m  Average  Wave  Height  Along  the 

Pier 

Direction 

(deg) 

Period 

(sec) 

Transmission  Coefficient  (HyH„) 

Maximum  Wave  Height,  m  (ft) 

Existing 

Dogleg 

Existing 

Dogleg 

221 

8 

0.50 

0.32 

3.0  (10.2) 

4.7  (15.4) 

10 

0.66 

0.45 

2.3  (7.5) 

3.4  (11.2) 

12 

0.79 

0.56 

1.9  (6.2) 

2.7  (8.9) 

14 

0.87 

0.67 

1.7  (5.6) 

2.3  (7.5) 

16 

0.92 

0.65 

1.7  (5.6) 

2.4  (7.9) 

20 

0.87 

0.84 

1.8  (5.9) 

1.8  (5.9) 

251 

8 

0.45 

0.25 

3.4  (11.2) 

6.1  (20.0) 

10 

0.58 

0.33 

2.6  (8.5) 

4.6  (15.1) 

12 

0.75 

0.42 

2.0  (6.6) 

3.6(11.8) 

14 

0.86 

0.53 

1.8  (5.9) 

2.9  (9.5) 

16 

0.90 

0.55 

1.7  (5.6) 

2.8  (9.2) 

20 

1.01 

0.64 

1.5  (4.9) 

2.4  (7.9) 

291 

8 

0.33 

0.21 

4.6  (15.1) 

7.2  (23.6) 

10 

0.48 

0.23 

3.1  (10.2) 

6.6  (11.8) 

12 

0.51 

0.29 

3.0  (9.8) 

5.2  (17.1) 

14 

0.54 

0.39 

2.8  (9.2) 

3.9  (12.8) 

Table  29 

Example  of  Average  Transmission  Coefficient  Calculation  for 

8-sec  Wave  from  221  deg 

Target 

Hmo  (m) 

Transmission  Coefficient  (H,/H„)  for  Gauge  Number 

1 

2 

3 

4 

Average 

2.5 

0.39 

0.19 

0.35 

0.38 

0.32 

3.5 

0.38 

0.17 

0.33 

0.40 

0.32 

4.5 

0.38 

0.17 

0.32 

0.41 

0.32 

Average 

0.38 

0.17 

0.33 

0.40 

0.32 

are  different  for  the  different  wave  heights,  the  results  are  close  in  magnitude  and 
show  that  averaging  the  results  from  different  incident  wave  heights  is  a  valid 
scheme  to  summarize  the  data. 
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Figure  53.  Importance  of  wave  height  on  wave  transmission  coefficient  for  a 
single  wave  period  and  direction  for  Hmo  wave  heights  of  2.5,  3.5,  and  4.5  m 


Another  method  was  used  to  further  compare  results  between  the  dogleg 
breakwater  and  existing  conditions.  This  technique  was  chosen  because  it  elimi¬ 
nates  small  differences  between  incident  deepwater  wave  conditions  when  waves  are 
generated  for  different  plans  or  from  different  directions.  Ideally  the  deepwater 
wave  conditions  would  always  be  identical;  however,  this  was  not  the  case.  The 
difference  between  transmission  coefficients  of  the  dogleg  breakwater  and  existing 
conditions  then  is  normalized  by  the  transmission  coefficients  for  existing  condi¬ 
tions.  The  formula  for  the  normalized  transmission  coefficients  is  given  as  follows: 


<  K  H, 

( _ l)  -( _ L) 

V  rr  ' doeles  V  T  T  ' 


K= 1  + 


H 


existing 


(— ) 

'  ' existing 


(16) 


in  which  Hi  represents  the  wave  height  of  the  ith  wave  gauge  and  H0  is  the  wave 
height  at  -29  m  (-95  ft).  Results  from  this  expression  show  the  increased  or 
decreased  effectiveness  of  the  dogleg  breakwater  over  existing  conditions.  The 
dogleg  breakwater  is  more  effective  in  reducing  the  wave  height  at  the  pier  than 
existing  conditions  if  the  result  is  less  than  unity.  If  the  results  are  unity,  there  is  no 
change  in  wave  height,  and  if  the  result  is  greater  than  unity,  the  dogleg  breakwater 
is  less  effective  than  existing  conditions.  The  decimal  percent  that  Na  is  above  or 
below  unity  is  how  much  less  or  more,  respectively,  the  dogleg  reduces  the  wave  at 
the  pier  than  the  existing  breakwater.  For  example,  given  the  same  deepwater  wave 
height  for  the  existing  condition  and  the  dogleg  breakwater  if  Na  =  0.75,  then  the 
wave  height  at  the  pier  with  the  dogleg  in  place  is  only  75  percent  as  high  as  it 
would  be  for  existing  conditions. 
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The  dogleg  breakwater  does  provide  significant  improvements  over  existing 
conditions.  The  dogleg  breakwater  and  existing  conditions  for  waves  from  221, 
251,  and  291  deg  referenced  to  north  are  in  Figure  54.  This  gives  an  indication  of 
the  decrease  in  wave  height  at  the  pier  with  the  dogleg  breakwater.  As  was  the  case 
with  the  transmission  coefficient  results,  the  normalized  transmission  coefficients 
represent  averaged  wave  heights  for  all  wave  cases  from  a  particular  direction. 

The  dogleg  breakwater  has  little  effect  on  wave  energy  reaching  gauge  2  for  waves 
from  291  deg.  Gauge  2  is  sheltered  from  waves  from  this  direction  and  increasing 
the  breakwater  would  have  little  effect. 


:§  Wave  Gage  Number 


a  Waves  from  221  deg 
Waves  from  251  deg 
-A-  Waves  from  291  deg 


Figure  54.  Normalized  transmission  coefficient  for  the  dogleg  breakwater 

The  straight  and  dogleg  breakwaters  were  compared  using  the  normalized 
transmission  coefficient  given  in  equation  16,  and  the  results  for  gauges  1-4  are 
shown  in  Table  30.  In  general,  the  dogleg  breakwater  provides  a  normalized  trans¬ 
mission  coefficient  of  0.65  while  the  straight  breakwater  provides  a  normalized 
transmission  coefficient  of  0.69.  Thus,  the  dogleg  breakwater  provides  slightly  more 
sheltering  than  the  straight  breakwater. 

Increased  Harbor  Usability  with  the  Dogleg 
Extension 


Before  considering  the  increased  harbor  usability  when  the  dogleg  breakwater  is 
present,  some  limitation  of  the  WIS  data  set  must  be  addressed.  Kaumalapau 
Harbor  is  directly  exposed  to  waves  propagating  from  180  deg  (due  south)  clock¬ 
wise  to  approximately  330  deg  (30  deg  west  of  north).  Therefore,  the  harbor  is 
exposed  to  waves  from  the  south  and  west  of  Lanai.  The  WIS  Pacific  Phase  I  grid  is 
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Table  30 

Comparison  Between  the  Straight  and  Dogleg  Breakwater 
Configurations  Using  the  Normalized  Transmission  Coefficients 
(Na) 

Direction 

Configuration 

Gauge  Number 

2 

3 

1 

4 

221 

Straight 

0.64 

0.76 

0.81 

0.77 

Dogleg 

0.62 

0.74 

0.80 

0.76 

251 

Straight 

0.55 

0.57 

0.65 

0.74 

Dogleg 

0.50 

0.51 

0.58 

0.70 

Average  Straight 

0.59 

0.67 

0.73 

0.76 

Average  dogleg 

0.56 

0.63 

0.69 

0.73 

Average  Straight  for  all  gauges  =  0.69. 

Average  dogleg  for  all  gauges  =  0.65. 

also  limited  in  size  to  the  west  and  south  of  Lanai,  as  shown  in  Figure  3.  The  WIS 
model  does  not  include  waves  that  are  generated  in  the  southern  hemisphere,  which 
are  waves  that  affect  Kaumalapau  Harbor  during  the  summer  months.  The  grid  size 
also  is  large,  and  low,  locally  generated  waves  are  not  included  in  the  wave  field 
calculated.  Limitations  of  the  WIS  Pacific  Phase  I  data  set  also  do  not  include 
effects  of  wave  sheltering  from  the  Hawaiian  Islands.  Additionally,  the  two  WIS 
stations  are  located  some  distance  from  Kaumalapau  Harbor  and  diffraction  of 
waves  around  Lanai  is  not  known. 

Given  the  limitations  of  the  WIS  data,  a  frequency  of  occurrence  table  was 
developed  (Table  31).  The  table  shows  frequency  of  occurrence  of  waves  over  a 
20-  year  record  from  January  1,  1956  through  January  1,  1976  from  wave  directions 
that  affect  Kaumalapau  Harbor.  Percent  of  occurrence  values  were  averaged  over 
the  two  WIS  stations.  The  last  four  columns  in  the  table  were  calculated  using 
transmission  coefficient  values  from  Table  28  and  the  calculated  deepwater  wave 
height,  period,  and  direction  from  the  20-year  WIS  record.  For  example,  waves 
from  206-235  occurred  only  0.05  percent  of  the  time.  If  existing  conditions  are 
present,  the  waves  along  the  wharf  are  reduced  to  less  than  1.5  m  (4.9  ft)  0.03  per¬ 
cent  of  the  20  years.  Since  the  dogleg  breakwater  provides  better  protection,  the 
waves  along  the  wharf  are  reduced  to  less  than  1.5  m  (4.9  ft)  0.04  percent  for  the 
20  years.  The  last  two  columns  indicate  the  percentage  of  time  wave  heights  along 
the  wharf  are  less  than  1.5  m  (4.9  ft)  for  waves  from  a  particular  direction  for  the 
existing  conditions  and  the  dogleg  extension.  Waves  from  291  deg  occur  4.8  per¬ 
cent  of  the  time  over  the  20-year  WIS  record.  For  2.09  percent  of  the  time,  the 
waves  will  be  from  291  deg  and  less  than  1.5  m  (4.9  ft)  along  the  wharf  for  the 
existing  conditions  while  for  the  dogleg  extension,  waves  meet  the  criteria 
3.71  percent  of  the  time.  Waves  meet  the  criteria  from  this  direction  44  percent  of 
the  time  for  the  existing  conditions  and  77  percent  of  the  time  for  the  dogleg 
extension. 
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Table  31 

Frequency  of  Occurrence  for  Waves  Used  in  Physical  Model 

Percent  of  Occurrence  Over  20-Year 
Period 

Percentage  of  Time 
Waves  from  a  Particular 
Direction  Generate 
Waves  at  the  Wharf 
Less  Than 

1.5  m  (4.9  ft) 

Wave 

Direction 

Bins 

Rep. 

Wave 

Direction 

WIS 

Data 

Meets  Wave 
Criteria  Along 
the  Pier  with 
Existing 
Conditions 

Meets  Wave 
Criteria 
Along  the 
Pier  with 
Dogleg 
Extension 

Existing 

Conditions 

Dogleg 

Extension 

206  -  235 

221 

0.05 

0.03 

0.04 

60 

80 

235  -  264 

251 

0.10 

0.08 

0.09 

80 

90 

265-315 

291 

4.8 

2.09 

3.71 

44 

77 

On  a  wave-by-wave  basis,  the  results  are  shown  in  Table  32.  The  numbers  in 
the  table  represent  the  combined  number  of  occurrences  of  a  3-hr  event  over  the 
20  years  of  WIS  data  for  the  two  stations.  The  unshaded  area  represents  maximum 
waves  that  will  meet  the  design  criteria  for  existing  conditions.  The  lightly  shaded 
area  represents  maximum  waves  that  will  meet  the  criteria  for  the  dogleg  extension, 
whereas  the  darkly  shaded  area  represents  waves  that  do  not  satisfy  the  criteria  of  a 
1.5-m  (4.9-ft)  wave  height  or  less  along  the  pier  face  with  the  dogleg  extension. 
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Table  32 

Number  of  Wave  Occurrences  from  20-year  WIS  Data  Set  as  a  Function  of  Wave  Parameters  and  the  Design 
Acceptance  Criteria  for  Waves  Along  the  Wharf 
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[Continued] 


10  Three-Dimensional  Stability 
Study1 


Armor  units  placed  on  the  ends  of  breakwaters  often  are  vulnerable  to  waves 
which  overtop  and  diffract  around  the  structure  because  the  units  are  subjected  to 
forces  from  the  opposing  direction  of  a  two-dimensional  situation.  Heavier  units 
often  are  required  for  breakwater  stability  on  the  head  of  the  structure.  To  deter¬ 
mine  stability  of  the  entire  Kaumalapau  breakwater,  three-dimensional  experiments 
were  conducted. 

The  Model 


Design  of  model 

Three-dimensional  stability  experiments  were  conducted  at  a  geometrically 
undistorted  linear  scale  of  1 :49,  model  to  prototype,  on  the  proposed  breakwater  at 
Kaumalapau  Harbor,  Hawaii.  Experiments  were  performed  on  the  same  model  used 
for  the  physical  harbor  model  experiments.  Scale  was  based  on  size  availability  of 
model  armor  units,  the  capabilities  of  the  available  wave  generator  to  produce 
required  wave  heights  at  the  modeled  water  depth,  and  necessary  size  to  also 
conduct  three-dimensional  harbor  model  experiments. 

Because  the  specific  weights  of  water  and  armor  layer  material  differed  between 
the  model  and  prototype,  the  transference  equation  of  Hudson  (1975)  was  used  to 
determine  the  model  scale  that  most  closely  represented  prototype  weights  for  the 
core-locs  available  at  WES: 


go, 

(K)P 


(ya)m 

(O 

3 

'(%  -  1' 

(Y  a)p 

UJ 

(SX  -  \ 

(17) 


in  which  the  subscripts  m  and  p  refer  to  model  and  prototype  quantities,  respec¬ 
tively,  and  ljlp  is  the  linear  scale  of  the  model. 


‘Written  by  Ernest  R.  Smith. 
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A  1:49  scale,  (ya)p  =  2,320  kg/m3  (145  lb/ft3),  and  model  core-loc  weights  of 
150  and  220  g  used  in  Equation  17  yielded  prototype  weights  of  18.1  and 
31.4  tonnes  (20  and  34.6  tons),  respectively,  and  were  used  as  primary  armor  layer 
Wj  (Figures  55  and  56).  Time  relations  were  scaled  according  to  Froude  Model 
Law  (Stevens  et.  al.  1942),  and  model-to-prototype  relations  were  derived  in  terms 
of  /  and  t  shown  in  Table  15. 


Figure  55.  Cross  section  of  Station  0  +  00 


Figure  56.  Cross  section  of  Station  2  +  00 
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Scale  effects  of  viscous  forces  associated  with  flow  through  the  underlayer  and 
core  of  the  proposed  breakwater  were  addressed  using  the  method  of  Keulegan 
(1973)  to  assure  that  flow  through  the  model  structure  was  turbulent.  Stone  weights 
of  18.6  to  29.0  g,  0.68  to  2.3  g,  and  29.0  to  60.0  g  with  a  specific  gravity  of 
approximately  2.68  represented  prototype  sizes  of  3.1  to  4.8  tonnes  (3.4  to 
5.3  tons),  0.5  to  68  kg  (1  to  150  lb),  and  4.8  to  10.0  tonnes  (5.3  to  1 1.0  tons)  for 
underlayer  stone  W2,  core  material  W3,  and  armor  cap  stone,  W4,  respectively. 

Model  underlayer,  core  material,  and  armor  cap  stone  remained  the  same  for  all 
three-dimensional  stability  plans,  and  are  listed  in  Figures  55  and  56  with  corre¬ 
sponding  prototype  sizes. 

Model  construction 

The  proposed  breakwater  was  constructed  on  the  existing  structure  at  a  crest 
elevation  of  +6. 1  m  (+20  ft)  mllw  and  a  crown  width  of  6. 1  m  (20  ft).  The  new 
breakwater  consisted  of  an  armor  layer,  an  underlayer,  and,  where  necessary,  a  core. 
Typical  cross  sections  of  the  breakwater  are  shown  in  Figures  55  and  56  at 
Stations  0+00  and  2+00,  respectively.  It  was  necessary  in  certain  locations  to 
excavate  material  from  the  existing  breakwater  to  place  the  armor  layer  and 
underlayer.  This  excavated  material  was  used  as  a  base  near  the  end  of  the  pro¬ 
posed  structure.  Core  material  was  added  only  in  locations  where  it  was  required  to 
raise  the  elevation  of  the  existing  breakwater  to  the  bottom  elevation  of  the 
underlayer. 

Construction  of  the  model  breakwater  simulated  prototype  construction  as 
closely  as  possible.  The  core  and  underlayer  material  were  dumped  by  shovel, 
smoothed  to  grade,  and  compacted  with  hand  trowels  to  simulate  consolidation  that 
would  have  occurred  due  to  wave  action. 

The  armor  layer  was  comprised  of  either  18. 1 -tonne  or  3 1.4-tonne  (20-  or 
34.6-ton)  core-locs  placed  on  a  IV  to  1.5H  slope.  Core-Iocs  extended  to  -10.7  m 
(-35  ft)  mllw  on  the  sea  side  of  the  structure,  and  extended  into  the  existing  break¬ 
water  profile  on  the  harbor  side  one  armor  unit  height.  Core-loc  armor  units  were 
placed  in  a  single  armor  layer  using  a  selective  random  placement  described  by 
Melby  and  Turk  (1995).  The  first  row  of  core-locs  were  aligned  with  vertical  flukes 
abutting  the  adjacent  units  (Figure  57).  The  second-row  units  were  placed  in  a 
manner  that  the  flukes  overlapped  the  waist  portion  of  the  first-row  units  (Fig¬ 
ure  57).  Units  above  these  were  placed  in  a  random  fashion  with  the  exception  that 
no  unit  was  placed  on  the  slope  with  vertical  flukes  directly  above  a  unit  also  placed 
with  vertical  flukes  (Figure  58).  Basic  two-dimensional  research  experiments  with 
core-locs  have  shown  that  two  units  placed  atop  each  other  (both  with  flukes 
oriented  vertically)  do  not  interlock  well  and  can  create  a  weakness  in  the  armor 
layer. 

The  number  of  core-loc  units  placed  on  the  breakwater,  or  density  of  units,  was 
based  on  the  equation: 

^V1  (18) 
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Figure  57.  Placement  of  toe  core-locs 


Figure  58.  Placement  of  core-locs 
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where  N  is  the  number  of  units  in  a  given  area  A,  V  is  the  armor  unit  volume,  and 
(p  is  the  packing  density  coefficient,  which  is  dependent  upon  armor  layer  thickness 
and  armor  layer  porosity.  Armor  layer  thickness  is  equal  to  about  0.92  of  the 
respective  core-loc  leg  length  and  the  average  armor  layer  porosity  is  about 
60  percent. 

The  core-loc  is  a  relatively  new  armor  unit  and  studies  are  ongoing  to  determine 
placement  density  of  the  units.  Results  of  basic  research  experiments  indicated  that 
the  packing  density  coefficient  for  core-locs  ranges  from  0.58  <<p<  0.63.  The 
lower  value,  <p  =  0.58,  (Plans  4-4E)  was  used  for  initial  buildings  of  the  model 
Kaumalapau  breakwater. 

A  rib  cap  was  placed  at  the  crown  of  the  structure  to  an  elevation  of  +6. 1  m 
(+20  ft)  mllw.  It  was  assumed  that  the  rib  cap  would  be  stable  in  the  prototype; 
therefore,  it  was  not  necessary  that  the  cap  be  dynamically  similar  to  the  prototype. 
The  model  rib  cap,  constructed  of  wood,  was  geometrically  similar  to  the  prototype 
and  was  anchored  in  the  model,  thus  ensuring  proper  transmission,  reflection,  and 
dissipation  of  wave  energy  and  the  assumed  stability  of  the  structure.  Individual 
ribs  were  0.9-m  (3-ft)  wide,  1.5-m  (5-ft)  high,  6.1-m  (20-ft)  long,  and  spaced  1.8  m 
(6  ft)  on  centers.  The  ribs  were  oriented  at  a  90-deg  angle  to  the  longitudinal  axis  of 
the  breakwater.  The  rib  cap  included  0.6-m-  (2-ft-)  wide  runners  placed  1.5  m  (5  ft) 
on  center  from  the  rib  ends  (Figure  59).  POD  desired  to  use  a  rounded  rib  cap  end 
at  the  breakwater  head;  however,  previous  experiments  with  rib  caps  indicated  that  a 
rounded  head  in  the  model  did  not  allow  for  interlocking  between  the  cap  and  armor 
units  (Bottin,  Markle,  and  Mize  1987).  Bottin,  Markle,  and  Mize  found  that  inter¬ 
locking  between  units  in  the  model  could  be  achieved  using  a  square-end  rib  cap. 
Therefore,  a  square-end  rib  cap  was  used  in  the  present  study,  but  a  rounded  rib  cap 
can  be  constructed  in  the  prototype  with  no  effect  on  stability  anticipated. 


Figure  59.  Dimensions  of  rib  cap 
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Study  facilities  and  equipment 

The  basin  used  for  physical  model  harbor  response  experiments  also  was  used  to 
conduct  the  stability  study.  Stability  experiments  were  performed  from  two  wave 
directions;  221  and  25 1  deg  referenced  to  north  which  were  selected  based  on  the 
direction  of  the  highest  waves  approaching  the  harbor.  The  head  of  the  proposed 
breakwater  is  exposed  to  waves  approaching  from  the  221 -deg  wave  direction,  and 
would  be  most  vulnerable  to  wave  attack  from  this  direction.  Therefore,  experi¬ 
ments  to  optimize  the  breakwater  were  initiated  from  the  221 -deg  direction.  After  a 
stable  plan  was  developed,  the  wave  generators  were  moved  to  the  25 1 -deg  wave 
direction  and  stability  was  verified  for  the  optimal  plan. 

Wave  conditions  studied  from  the  221-  and  25 1 -deg  directions  are  given  in 
Table  33,  which  lists  the  wave  periods  used  during  the  stability  study,  zero-moment 
deepwater  wave  height  Ha,  zero-moment  wave  height  at  the  wave  generator  depth 
of  29  m  (95  ft)  mllw  H\  average  zero-moment  wave  height  at  the  structure  Hd,  and 
duration  of  each  wave  series  in  model  seconds.  The  design  wave  was  specified  in  a 
reconnaissance  report  as  a  9.8-sec,  8.5-m  (28-ft)  deepwater  wave  (USAED, 
Honolulu  1993).  To  determine  a  factor  of  safety  for  the  structure,  waves  exceeding 
the  design  height  were  generated  for  9.8-  and  12-sec  periods.  Wave  heights 
associated  with  16-sec  periods  from  the  221-  and  251 -deg  directions  are  low  in 
magnitude;  therefore,  waves  up  to  Ha=  8.5  m  (28  ft)  were  generated  for  16-sec 
waves.  The  total  duration  of  the  conditions  shown  in  Table  33  represents  a 
prototype  storm  of  approximately  24  hr.  All  wave  conditions  were  conducted  at  a 
storm  surge  water  level  of +1.5  m  (+5.0  ft)  mllw.  Because  the  breakwater  location 
is  in  relatively  deep  water,  a  breaking  wave  condition  was  not  reached  during 
calibration  without  the  structure  in  place.  However,  because  of  the  interaction  of 
reflected  and  incident  waves,  breaking  waves  were  observed  during  stability 
experiments  with  the  structure  in  place. 

Study  procedures 

Photographs  were  taken  of  the  harbor  side,  head,  and  sea  side  before  the 
experiment  was  initiated  without  water  in  the  basin.  Following  before-experiment 
photographs,  the  basin  was  flooded  to  +1.5  m  mllw  and  tire  structure  was  exposed 
to  9.8-sec,  3.0-m  (10-ft)  waves.  The  low-level  series  allowed  settling  and  nesting  of 
the  newly  constructed  section  which  would  occur  under  typical  daily  wave 
conditions  prior  to  being  exposed  to  a  design  level  storm.  The  remainder  of  the 
storm  conditions  for  the  particular  direction  were  generated  upon  completion  of  the 
low-level  waves.  Response  of  the  structure  was  recorded  during  and  after  each 
wave  condition.  A  detailed  inspection  of  the  structure  was  performed  between  wave 
series,  and  effects  of  the  waves  on  individual  units,  toe  stability,  and  the  general 
condition  of  the  breakwater  were  recorded.  The  basin  was  drained  upon  completion 
of  the  entire  storm  series  or  after  the  structure  had  been  severely  damaged,  and  after 
test  photographs  were  taken  of  the  harbor  side,  head,  and  sea  side.  The  same 
procedure  was  followed  if  the  plan  was  subjected  to  additional  storm  conditions. 
Before  and  after  photographs  are  located  in  Appendix  B. 
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Table  33 

Kaumalapau  Stability  Study  Conditions 

TP 

sec 

H„ 

m  (ft) 

H’ 

m  (ft) 

H„  (221  deg) 
m(ft) 

H„  (251  deg) 
m  (ft) 

Duration 

sec 

(model) 

9.8 

3.6  (11.8) 

3.0  (10.0) 

3.2  (10.5) 

3.0  (10.0) 

880 

9.8 

5.5  (18.0) 

4.6  15.0) 

4.8  (15.7) 

4.4  (14.4) 

880 

9.8 

7.3  (24.0) 

6.1  (20.0) 

6.4  (21.0) 

5.9  (19.4) 

880 

9.8 

8.51  (28.0) 

7.2  (23.6) 

7.5’(24.7) 

6.9'  (22.6) 

880 

9.8 

9.1  (30.0) 

7.6  (25.0) 

8.3  (27.3) 

7.3  (24.0) 

880 

9.8 

9.9  (32.5) 

8.3  (27.2) 

8.7  (28.5) 

8.0  (26.2) 

880 

12.0 

5.5  (18.0) 

4.6(15.0) 

4.3  (14.1) 

4.3  (14.1) 

880 

12.0 

7.3  (24.0) 

6.1  (20.0) 

5.8  (19.0) 

5.7(18.7) 

880 

12.0 

8.5  (28.0) 

7.2  (23.6) 

6.8  (22.3) 

6.7  (22.0) 

880 

12.0 

9.1  (30.0) 

7.6  (25.0) 

7.3  (24.7) 

7.2  (23.6) 

880 

12.0 

10.8(35.4) 

9.1  (30.0) 

8.6  (28.2) 

8.5  (28.0) 

880 

16.0 

4.9  (16.0) 

4.6  (15.0) 

3.9  (12.8) 

4.5  (14.8) 

880 

16.0 

6.6(21.7) 

6.1  (20.0) 

5.2  (17.1) 

6.0  (19.7) 

880 

16.0 

8.5  (28.0) 

7.9  (26.0) 

6.7  (22.0) 

7.7  (25.3) 

880 

Note: 

H0  -  deepwater  zero-moment  wave  height 

H’  -  zero-moment  wave  height  at  wave  generator  depth  (95  ft  prototype) 

Hd  -  average  zero-moment  wave  height  at  structure 

1  Design  wave  height  determined  by  Pacific  Ocean  Division. 

Visual  inspections  were  made  during  and  after  wave  action  on  the  structure. 
Qualitative  results  were  taken  for  overtopping,  and  the  wave  height  and  category  of 
overtopping  (minor,  moderate,  or  major)  were  noted  for  each  wave  condition. 

Results 

Results  from  the  harbor  model  experiments  indicated  that  the  dogleg  breakwater 
configuration,  Plan  4,  was  the  more  desirable  plan;  therefore,  all  stability  experi¬ 
ments  were  conducted  with  the  dogleg  configuration. 

Nine  different  stability  plans  were  studied  which  differed  by  armor  unit  weight, 
armor  placement  density,  and  toe  protection  schemes.  The  following  paragraphs 
summarize  each  plan  studied.  Unless  otherwise  specified,  wave  heights  are  referred 
to  as  wave  height  at  the  structure  Hd 
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Plan  4 


Plan  4  armor  consisted  of  all  18. 1-tonne  (20-ton)  core-locs  (Photos  B1  through 
B3).  Experiments  were  initiated  with  a  low-level  (9.8-sec,  3.2-m  (10.5-ft))  wave 
spectra  from  the  221-deg  direction  for  a  model  duration  of  880  sec  to  allow  settling 
and  consolidation  of  the  armor  and  underlayer  materials.  Following  the  low-level 
waves,  heights  ( Hd )  of  4.8,  6.4,  7.5,  and  8.3  m  (15.8, 21.0,  24.7,  and  27.3  ft, 
respectively)  were  generated  at  9.8  sec.  The  breakwater  was  stable  for  4.8-m 
(15.8  -ft)  waves,  but  six  units  were  displaced  off  the  sea  side  of  the  structure  at  the 
shore  interface  during  6.4-m  (2 1.0-ft)  waves.  The  structure  continued  to  unravel  in 
this  area  for  the  design  height  of  7.5  m  (24.7  ft)  (20  additional  units  displaced),  and 
for  the  subsequent  series  of  8.3-m  (27.2-ft)  waves  (28  total  units  displaced. 

Photo  B4).  However,  the  remainder  of  the  breakwater  was  stable  with  no  units  dis¬ 
placed  for  the  conditions  studied  (Photos  B5  and  B6).  The  damaged  area  is  located 
in  shallow  water  on  a  reef.  The  alignment  of  the  breakwater  with  the  underwater 
bathymetry  of  the  reef  causes  wave  energy  to  converge  at  the  base  of  the  structure, 
and  incident  waves  from  deeper  water  break  in  this  region. 

Plan  4A 

The  damaged  area  was  repaired  and  additional  18.1-tonne  (20-ton)  core-locs 
were  placed  on  the  sea  side  at  the  shoreward  end  of  the  structure.  The  1 8. 1-tonne 
(20-ton)  core-locs  were  placed  farther  shoreward  to  the  seawall  in  an  attempt  to  pro¬ 
vide  additional  support.  Experiments  with  Plan  4A  were  initiated  with  9.8-sec, 

4.8- m  (15.8 -ft)  waves.  Six  units  were  displaced  on  the  sea  side  at  the  shoreward 
end,  and  subsequent  waves  of  6.4  m  (21  ft)  displaced  20  additional  units 

(Photo  B7).  The  remainder  of  the  structure  was  stable  during  Plan  4A  experiments. 

Plan  4B 

To  provide  additional  armor  in  the  problem  area,  18.1 -tonne  (20-ton)  core-locs 
were  placed  farther  offshore  for  Plan  4B  to  provide  a  wider  armor  area  and  gentler 
slope  (Photo  B8).  Experiments  of  Plan  4B  were  initiated  with  9.8-sec,  4.8-m 
(15.8-ft)  waves,  which  displaced  six  units.  The  wider  toe  area  was  not  sufficient  to 
stabilize  the  shoreward  end,  and  was  damaged  after  9.8-sec,  6.4-m  (21.0-fit)  waves 
(Photo  B9).  It  appeared  from  Plan  4B  experiments  that  the  toe  was  in  a  region  of 
high  energy,  and  units  above  the  toe  were  not  displaced  directly  by  wave  energy,  but 
by  toe  failure. 

Plan  4C 

To  avoid  placing  the  toe  units  in  the  region  of  high  energy,  a  narrower  armor 
section  of  18. 1-tonne  (20-ton)  core-locs  was  constructed  at  the  shoreward  end  of  the 
sea  side  for  Plan  4C  (Photo  BIO).  However,  the  plan  also  failed  at  the  toe  after 

9.8- sec,  6.4-m  (21.0-ft)  waves  were  generated  (Photo  B1 1). 
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Plan  4D 


To  protect  the  toe  units  on  the  sea  side  at  the  shoreward  end,  a  15.2-m-long 
(50-ft-long)  spur  of  18. 1 -tonne  (20-ton)  core-locs  was  constructed  perpendicular  to 
the  breakwater.  Plan  4D  (Photo  B 12).  The  spur  was  destroyed,  and  the  toe  of 
Plan  4D  failed  after  9.8-sec,  6.4-m  (21.0-lit)  waves  (Photo  B13).  Experiments  with 
Plans  1  through  ID  showed  that  the  area  of  toe  instability  on  the  sea  side  at  the 
shoreward  end  was  subjected  to  high  wave  forces  and  velocities.  The  area  is  located 
in  shallow  water,  and  the  breakwater  height  is  low  relative  to  the  depth.  Toe  units  of 
the  breakwater  in  this  area  do  not  have  sufficient  weight  to  support  unraveling  of  the 
structure  toe. 

Plan  4E 

To  stabilize  the  toe  in  the  problem  area,  a  strip  of  sheet  metal  was  placed  in  front 
of  the  toe  units  on  the  reef,  Plan  4E,  to  simulate  a  toe  trench  in  prototype 
(Photo  B 14).  The  sheet  metal  toe  strip  originated  at  the  shoreward  end  of  the 
structure  and  extended  down  the  reef  for  an  overall  length  of  41. 1  m  (135  ft), 
prototype.  Plan  4E  was  subjected  to  the  series  of  9.8-sec  waves  listed  in  Table  33. 
The  18. 1 -tonne  (20-ton)  core-locs  in  the  vicinity  of  the  toe  trench  were  stable; 
however,  three  units  were  displaced  off  the  breakwater  head  during  7.5-m  (24.7-ft) 
waves.  The  breakwater  was  subjected  to  12-sec,  4.3-  and  5.8-m  (14.2-  and  19.0-ft) 
waves  and  16-sec,  3.9-  and  5.2-m  (12.8-  and  17.0-ft)  waves.  No  additional  core- 
locs  were  displaced  off  the  structure  for  12-sec  waves  or  the  16-sec,  3.9-m  (12.8-ft) 
condition;  however,  units  were  observed  to  move  on  the  head  during  wave  action. 

The  sea  side  of  the  structure  was  damaged  during  16-sec,  5.2-m  (17.0-ft)  waves. 

Plan  4E  was  rebuilt  with  1 8. 1-tonne  (20-ton)  core-locs  (Photos  B15  through 
B17),  and  experiments  were  conducted  for  9.8-sec  waves  up  to  7.5  m  (24.7  ft).  The 
overall  condition  of  the  breakwater  was  good,  but  the  head  was  damaged;  four  core- 
locs  were  displaced  off  the  head  near  the  crown  and  the  underlayer  and  stone 
beneath  the  rib  cap  was  exposed  (Photos  B18  through  B20).  Experiments  of 
Plan  4E  continued  with  12-sec  waves  up  to  6.8  m  (22.3  ft).  The  head  deteriorated 
further,  but  no  additional  displacement  of  units  occurred  and  the  breakwater  con¬ 
dition  was  good  (Photos  B21  through  B23).  The  plan  was  subjected  to  16-sec  waves 
for  all  heights  shown  in  Table  10-1.  Four  additional  units  were  displaced  on  the  sea 
side  of  the  structure  during  5.2-m  (17. 1-ft)  waves,  and  five  additional  units  were 
displaced  on  the  sea  side  during  6.7-m  (22.0-ft)  waves  (Photos  B24  through  B26). 
Plan  4E  was  then  subjected  to  9.8-sec,  8.7-m  (28.7-ft)  waves,  and  the  head  was 
severely  damaged  (Photos  B27  through  B29). 

Plan  4F 

The  breakwater  head  was  rebuilt  using  18.1 -tonne  (20-ton)  core-locs  and  a 
denser  placement,  (p=  0.63,  and  experiments  were  conductedfor  all  9.8-sec  waves. 
Units  on  the  head  began  to  settle  during  the  design  wave,  but  no  core-locs  were 
displaced  off  the  structure.  One  core-loc  was  displaced  during  8.7-m  (28.5-ft) 
waves.  The  entire  structure  was  rebuilt  using  the  denser  placement  of  units 
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(Photos  B30  through  B32).  Stability  experiments  of  Plan  4F  were  conducted  for  all 
periods  up  to  H0  =  8.5  m  (28.0  ft).  Two  units  were  displaced  for  9.8-sec  waves  at 
the  design  wave  height,  Hd  =  7.5  m  (24.7  ft)  (Photos  B33  through  B35).  No  further 
displacement  occurred  during  12-sec  waves,  but  units  were  observed  to  move  on  the 
head  (Photos  B36  through  B38).  One  additional  core-loc  was  displaced  during 
16-sec,  3.9-m  (12.8-ft)  waves,  but  no  displacement  occurred  for  5.2-m  (17.1-ft) 
waves.  Eleven  additional  core-locs  were  removed  by  16-sec,  6.7-m  (22.0-ft)  waves 
on  the  sea  side  of  the  breakwater  near  the  shoreward  end  (Photos  B39  through  B4 1). 
Rundown  of  waves  displaced  units  at  the  toe,  which  were  not  protected  by  the 
trench,  and  caused  instability  on  the  slope.  Waves  were  generated  for  9.8-  and 
12.0-sec  periods  for  heights  greater  than  the  design  height.  No  additional  units  were 
displaced  for  9.8-sec,  8.3-m  (27.3 -ft)  waves,  but  four  units  were  displaced  off  the 
head  for  12-sec,  7.3-m  (23.8-ft)  waves.  The  9.8-sec,  8.7-m  (28.5-ft)  wave  condition 
removed  six  additional  units  off  the  head,  and  the  12-sec,  8.6-m  (28.3-ft)  wave 
severely  damaged  the  harbor  side  of  the  structure  (Photos  B42  through  B44). 

Plan  4G 


Stability  experiments  of  the  breakwater 
continued  from  the  221 -deg  wave  direction 
with  Plan  4G  which  consisted  of  18. 1 -tonne 
(20-ton)  core-locs  placed  (<j)  =  0.62)  on  the  sea¬ 
side  and  harbor-side  trunks  and  3 1.4-tonne 
(34.6-ton)  units  placed  (<j)  =  0.62)  on  the  head 
of  the  structure  (Figure  60,  Photos  B45  through 
B47).  Plan  4G  was  subjected  to  the  conditions 
listed  in  Table  33  for  9.8-  and  12.0-sec  waves 
up  to  Ha  =  8.5  m  (28.0  ft),  and  16.0-sec  waves 
up  to  Ha  =  6.6  m  (21.6  ft).  Wave  hindcast  data 
indicated  that  waves  associated  with  the 
16.0-sec  period  from  the  221 -deg  direction  did 
not  approach  deepwater  heights  of  the  magni¬ 
tude  of  8.5  m  (28.0  ft);  therefore,  this  condition 
was  excluded  from  further  stability  experiments 
from  the  221 -deg  direction.  The  structure  was 
not  damaged  for  the  conditions  described  above 
(Photos  B48  through  B50).  Waves  were 
generated  for  heights  greater  than  H0  =  8.5  m 
(28  ft).  Two  31.4-tonne  (34.6-ton)  core-locs 
were  displaced  off  the  head  for  9.8-sec,  8.7-m 
(28.5-ft)  waves  (Photos  B5 1  and  B52).  The  harbor  side  of  the  breakwater  was 
severely  damaged  for  12-sec,  8.6-m  (28.2-ft)  waves  (Photos  B5 1  and  B52),  but  the 
sea  side  of  the  structure  remained  stable  for  all  waves  (Photo  B53). 


Figure  60.  Plan  4G 


Plan  4H 

To  determine  if  the  harbor  side  could  be  stabilized,  3 1.4-tonne  (34.6-ton)  core- 
locs  were  placed  (4>  =  0.62)  on  both  the  head  and  the  trunk  of  the  harbor  side 
(Figure  61,  Photos  B54  and  B55).  Due  to  the  lack  of  an  adequate  number  of 
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3 1.4-tonne  (34.6-ton)  model  core-iocs,  the 
breakwater  sea  side  was  constructed  (<j)  = 

0.62)  using  18.1-tonne  (20-ton)  core-locs 
(Photo  B56).  The  structure  was  subjected  to 
all  221 -deg  wave  conditions  listed  in 
Table  33  with  the  exception  of  16.0-sec, 

6.7-m  (22.0-ft)  waves.  Three  3 1.4-tonne 
(34.6-ton)  core-locs  were  displaced  off  the 
harbor  side  during  12-sec,  8.6-m  (28.2-ft) 
waves,  but  the  breakwater  was  in  stable 
condition  (Photos  B57  through  B59). 

The  wave  generator  was  moved  to  the 
251 -deg  wave  direction.  Plan  4H  was  rebuilt 
using  <j)  =  0.62  (Photos  B60  through  B62) 
and  subjected  to  9.8-  and  12.0-sec  waves  up 
to  the  design  wave  height  listed  in  Table  33, 

Hd  =  6.9  and  6.7  m  (22.6  and  22.0  ft) 
respectively,  and  16.0-sec  waves  up  to  Hd  = 

6.0  m  (19.7  ft).  One  18.1-tonne  (20-ton) 
core-loc  was  displaced  on  the  sea  side  near  Figure  61 .  Plan  4H 
the  toe  trench  during  9.8-sec,  6.9-m  (22.6-ft) 

waves,  but  the  breakwater  was  stable  (Photos  B63  through  B65).  Experiments 
continued  with  9.8-sec  waves  up  to  8.0  m  (26.2  ft)  and  12-sec  waves  up  to  8.5  m 
(28.0  ft)  (Table  33).  The  structure  was  stable  for  the  conditions  studied;  however, 
five  18.1 -tonne  (20-ton)  units  were  displaced  on  the  sea  side  and  one  3 1.4-tonne 
(34.6-ton)  core-loc  was  displaced  on  the  harbor  side  (Photos  B66  through  B68). 

The  breakwater  was  subjected  to  the  16.0-sec,  7.7-m  (25.3-ft)  wave  (Table  33) 
without  rebuilding,  and  no  further  damage  occurred  (Photos  B69  through  B71). 

Plan  4H  was  rebuilt  (cj)  =  0.62)  and  stability  experiments  were  conducted  for  9.8- 
and  12.0-sec  waves  up  to  H0  =  8.5  m  (28.0  ft),  and  16.0-sec  waves  up  to  Ha  =  6.6  m 
(21.7  ft),  and  no  damage  was  observed  (Photos  B72  through  B74).  The  structure 
remained  stable  when  subjected  to  waves  exceeding  the  design  height  for  9.8-sec 
and  12.0-sec  waves  (Table  33),  and  the  16.0-sec  wave  at  the  design  height 
(Photos  B75  through  B77). 

Summary 

Nine  stability  plans  were  studied  for  the  proposed  dogleg  breakwater.  It  was 
determined  that: 

a.  A  toe  trench  at  the  shoreward  end  on  the  sea  side  of  the  breakwater  was 

necessary  to  obtain  stability.  The  trench  had  an  overall  length  of  approxi¬ 
mately  41. 1  m  (135  ft),  prototype. 
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b.  Densely  placed  (^=0.63)  18. 1-tonne  (20-ton)  core-locs  were  stable  for  the 
design  wave  conditions,  but  failed  for  wave  heights  exceeding  the  design 
conditions. 

c.  A  heavier  core-loc,  3 1.4  tonnes  (34.6-tons),  placed  on  the  head  and  harbor 
side  at  (p=  0.62  was  found  to  be  stable  for  wave  heights  exceeding  the 
design  condition  from  the  22 1  -  and  25 1  -deg  wave  directions. 

Although  experiments  were  not  conducted  with  3 1.4-tonne  (34.6-ton)  units  on  the 
sea  side  of  the  breakwater  because  of  the  limited  number  of  model  units,  3 1.4-tonne 
(34.6-ton)  core-locs  placed  on  the  sea  side  are  expected  to  be  stable  for  the 
conditions  studied.  It  may  be  desirable  to  place  one  size  unit  on  the  entire  structure 
to  reduce  costs  of  constructing  an  additional  armor  unit  mold.  Also,  less  3 1.4-tonne 
(34.6-ton)  units  would  be  required  to  cover  an  area  than  18. 1-tonne  (20-ton)  units. 
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1 1  Conclusions  and 
Recommendations1 


Harbor  Response  to  Wind  Waves,  Swell,  and  Long 
Waves 

Field  data  were  collected  to  validate  physical  and  numerical  models  and  assess 
hindcast  data.  A  numerical  model  study  was  used  to  assist  in  optimizing  the 
placement  of  field  wave  gauges,  maximize  efficiency  of  physical  model  experi¬ 
ments  by  identifying  the  most  promising  modification  plans  and  most  damaging 
incident  wave  directions,  and  to  evaluate  harbor  oscillation  characteristics  of  the 
proposed  plan.  An  undistorted,  three-dimensional  physical  model  of  Kaumalapau 
Harbor  was  used  to  evaluate  harbor  response  to  short  period  waves  for  existing 
conditions  and  two  alternative  harbor  layout  configurations,  straight  and  dogleg 
breakwater  plans,  respectively,  for  three  incident  wave  directions  as  well  as 
optimize  the  stability  of  the  proposed  breakwater  design. 

Conclusions 

a.  For  existing  conditions,  the  area  just  west  of  the  west  end  of  the  wharf  is 
significantly  sheltered.  Wave  heights  in  this  area  can  be  expected  to  be 
less  than  half  the  incident  height  for  most  incident  conditions. 

b.  The  western  end  of  the  wharf  is  also  relatively  sheltered  if  compared  to 
the  middle  and  eastern  part  of  the  wharf  for  most  wave  periods  for 
existing  conditions. 

c.  The  longest  wave  (swell)  periods,  especially  those  from  the  northwesterly 
incident  directions,  are  amplified  at  the  middle  and  western  part  of  the 
wharf  but  not  at  the  eastern  part  for  existing  conditions.  This  result  is 
interesting  in  light  of  local  reports  that  high  northerly  swell  can  cause 
problems  in  Kaumalapau  Harbor. 
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d.  For  existing  conditions,  harbor  response  varies  considerably  with  incident 
wave  direction. 

e.  A  45.7-m  (150-ft)  straight  extension  to  the  existing  breakwater  will 
reduce  short-period  wind  wave  heights  at  the  wharf  to  about  60  to 
80  percent  of  those  in  the  existing  harbor. 

/.  The  45.7-m  (150-ft)  dogleg  breakwater  extension  is  slightly  more  effec¬ 
tive  than  the  straight  breakwater  for  reducing  short-period  wave  action  at 
the  pier  and  for  this  reason  is  the  recommended  plan. 

g.  Neither  the  straight  nor  the  dogleg  extension  met  short-period  wave 
reduction  criteria  established  by  POD.  Further  wave  reduction  could  be 
accomplished  but  the  breakwater  would  not  meet  the  navigational  criteria 
and  economic  constraints. 

h.  Although  the  dogleg  plan  is  prone  to  stronger  oscillations  than  the  existing 
harbor,  the  differences  will  have  little  operational  significance.  Wind 
waves  and  swell  can  be  expected  to  continue  as  a  primary  operational 
concern  which  will  probably  outweigh  any  practical  problems  due  to 
increased  long-period  oscillations. 

i.  Wave  absorbers  along  the  northeast  or  south  shores  of  the  embayment 
have  minimal  impact  on  wave  heights  at  the  pier. 

Nine  stability  plans  were  studied  using  the  dogleg  breakwater  configuration. 

Based  on  model  experiments,  it  was  determined  that: 

a.  A  toe  trench  at  the  shoreward  end  on  the  sea  side  of  the  breakwater  was 
necessary  to  obtain  armor  stability.  The  trench  had  an  overall  length  of 
approximately  41.1  m  (135  ft),  prototype. 

b.  Densely  placed  (0=0.63)  18. 1-tonne  (20-ton)  core-locs  were  stable  for 
design  wave  conditions,  but  failed  for  wave  heights  exceeding  the  design 
conditions. 

c.  A  heavier  core-loc,  31.4  tonnes  (34.6  tons),  placed  on  the  head  and 
harbor  side  (0=0.62),  was  found  to  be  stable  for  wave  heights  exceeding 
the  design  condition  from  the  221-  and  251 -deg  wave  directions. 

d.  It  may  be  desirable  to  place  one  size  unit  on  the  entire  structure  to  reduce 
costs  of  constructing  an  additional  armor  unit  form.  Additionally,  less 
armor  units  would  be  required  for  placement  if  the  breakwater  rehabilita¬ 
tion  was  constructed  entirely  of  the  larger  unit.  Therefore  it  is  recom¬ 
mended  that,  if  the  breakwater  is  constructed  using  one  unit,  31.4-tonne 
(34.6-ton)  core-locs  be  placed  using  0=0.62. 
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Figure  A1 .  Prototype  data  summary  for  January,  1 994,  outside  harbor 
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Figure  A2.  Prototype  data  summary  for  January,  1994,  inside  harbor 
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Figure  A5.  Prototype  data  summary  for  March,  1994,  outside  harbor 
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Figure  A7.  Prototype  data  summary  for  April,  1 994,  outside  harbor 
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Figure  A8.  Prototype  data  summary  for  April,  1994,  inside  harbor 
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Figure  A1 0.  Prototype  data  summary  for  May,  1 994,  inside  harbor 
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Figure  A1 1 .  Prototype  data  summary  for  June,  1 994,  outside  harbor 
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Figure  A14.  Prototype  data  summary  for  July,  1994,  inside  harbor 
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Figure  A15.  Prototype  data  summary  for  August,  1994,  outside  harbor 
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Photo  B2.  Plan  4,  breakwater  head  before  experiments 
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Photo  B4.  Sea-side  view  of  Plan  4  after  9.8-sec,  7.5-m  waves 
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Photo  B5.  Plan  4,  leeside  view  of  breakwater  after  9.8-sec,  7.5-m  waves 
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Photo  B6.  Head  view  of  Plan  4  after  9.8-sec,  7.5-m  waves 
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Photo  B9.  Sea-side  shoreward  end  view  of  Plan  4B,  after  9.8-sec,  6.4-m  waves 


Photo  BIO.  Sea-side  shoreward  end  view  of  Plan  4C  before  experiments 
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Photo  B1 2.  Sea-side  shoreward  end  view  of  Plan  4D  before  experiments 
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Photo  B13.  Sea-side  shoreward  end  view  of  Plan  4D,  after  9.8-sec,  6.4-m  waves 
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Photo  B16.  Plan  4E,  breakwater  head  before  experiments 
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Photo  B1 7.  Sea-side  view  of  Plan  4E,  before  experiments 
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Photo  B18.  Leeside  view  of  Plan  4E,  after  9.8-sec,  7.5-m  waves 
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Photo  B20.  Sea-side  view  of  Plan  4E  after  9.8-sec,  7.5-m  waves 
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Photo  B22.  Plan  4E,  breakwater  head  after  12-sec,  6.8-m  waves 
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Photo  B25.  Plan  4E,  breakwater  head  after  1 6-sec,  6.7-m  waves 
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Photo  B26.  Sea-side  view  of  Plan  4E  after  16-sec,  6.7-m  waves 
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Photo  B28.  Plan  4E,  breakwater  head  after  9,8-sec,  8.7-m  waves 
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Photo  B29.  Sea-side  view  of  Plan  4E  after  9.8-sec,  8.7-m  waves 
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Photo  B31 .  Plan  4F,  breakwater  head  before  experiments 


Photo  B32.  Sea-side  view  of  Plan  4F  before  experiments 
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Photo  B34,  Plan  4F,  breakwater  head  after  9.8-sec,  7.5-m  waves 
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Photo  B36.  Leeside  view  of  Plan  4F  after  12-sec,  6.8-m  waves 
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Photo  B37.  Plan  4F,  breakwater  head  after  12-sec,  6.8-m  waves 


Photo  B38.  Sea-side  view  of  Plan  4F,  after  1 2-sec,  6.8-m  waves 


B20 


Appendix  B  Three-Dimensional  Stability  Photographs 


kaumaUWau 

STABILITY 

plan  If  & 

,\FT£R 

i6s  a#w 

C7J5  ........... 


Photo  B40.  Plan  4,  breakwater  head  after  1 6-sec,  6.7-m  waves 
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Photo  B41 .  Sea-side  view  of  Plan  4F,  after  1 6-sec,  6.7-m  waves 
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Photo  B44.  Sea-side  view  of  Plan  4F  after  maximum  9.8-  and  12-sec  waves 
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Photo  B46.  Plan  4G,  breakwater  head  before  experiments 
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Photo  B50.  Sea-side  view  of  Plan  4G  after  9.8-  and  1 2-sec  waves  at  design  height  and  1 6-sec  waves 
up  to  H0  =  6.6  m 
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Photo  B53.  Sea-side  view  of  Plan  4G  after  maximum  9.8-  and  1 2-sec  waves 
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Photo  B57.  Leeside  view  of  Plan  4H  after  maximum  9.8-  and  1 2-sec  waves  and  1 6-sec  waves  up  to 
H0  =  6.6  m 
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Photo  B58.  Plan  4H,  breakwater  head  after  maximum  9.8-  and  1 2-sec  waves  and  1 6-sec  waves  up  to 
Ha  =  6.6  m 
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Photo  B62.  Sea-side  view  of  Plan  4H  before  experiments  from  251 -deg  direction 
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Photo  B63.  Leeside  view  of  Plan  4H  after  9.8-  and  1 2-sec  waves  at  design  height  and  1 6-sec  waves  up 
to  H0  =  6.6  m  from  251 -deg  direction 
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Photo  B64.  Plan  4H,  breakwater  head  after  9.8-  and  1 2-sec  waves  at  design  height  and  1 6-sec  waves 
up  to  H0  =  6.6  m  from  251 -deg  direction 
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Photo  B65. 


and  12-sec  waves  from  251 -deg  direction 
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Sea-side  view  of  Plan  4H  after  9.8-  and  12-sec  waves  at  design  height  and  16-sec  waves 
up  to  Ha  =  6.6  m  from  251  -deg  direction 


B34 


Leeside  view  of  Plan  4H  after  maximum  9.8- 
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Photo  B68.  Sea-side  view  of  Plan  4H  after  maximum  9.8-  and  1 2-sec  waves  from  251  -deg  direction 
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Photo  B69.  Leeside  view  of  Plan  4H  after  16-sec,  7.7-m  waves  from  251 -deg  direction 


Photo  B70.  Plan  4H,  breakwater  head  after  1 6-sec,  7.7-m  waves  from  251  -deg  direction 
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Photo  B71 .  Sea-side  view  of  Plan  4H  after  1 6-sec,  7.7-m  waves  from  251  -deg  direction 


Photo  B72.  Leeside  view  of  Plan  4H  after  9.8-  and  12-sec  waves  at  design  height  and  16-sec  waves  up 
to  Ha  =  6.6  m  from  251  -deg  direction,  repeat  test 
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Photo  B73.  Plan  4H,  breakwater  head  after  9.8-  and  1 2-sec  waves  at  design  height  and  1 6-sec  waves 
up  to  H0  =  6.6  m  from  251 -deg  direction,  repeat  test 


Photo  B74.  Sea-side  view  of  Plan  4H  after  9.8-  and  1 2-sec  waves  at  design  height  and  1 6-sec  waves 
up  to  H0  =  6.6  m  from  251  -deg  direction,  repeat  test 
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Photo  B75. 


Leeside  view  of  Plan  4H  after  maximum  9.8-  and  12-sec  waves  and  16-sec  waves  up  to 
H0  =  8.5  m  from  251 -deg  direction,  repeat  test 


Photo  B76.  Plan  4H,  breakwater  head  after  9.8-  and  1 2-sec  waves  and  1 6-sec  waves  up  to  H0  =  8.5  m 
from  251 -deg  direction,  repeat  test 
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Photo  B77.  Sea-side  view  of  Plan  4H  after  maximum  9.8-  and  1 2-sec  waves  and  1 6-sec  waves  up  to 
Ha  =  8.5  m  from  251 -deg  direction,  repeat  test 
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Notation 


a  Wave  amplitude,  m  (ft) 

o,  Incident  wave  amplitude,  m  (ft) 

Amp  Wave  amplification  factor 

(Aamp)'ff  Effective,  or  spectral,  wave  amplification  factor 

c  Wave  phase  speed,  m/sec  (ft/sec) 

cg  Wave  group  speed,  m/sec  (ft/sec) 

d  Water  depth,  m  (ft) 

dfar  Water  depth,  m  (ft) 

fp  Peak  spectral  frequency,  sec"1 

H  Wave  height,  m  (ft) 

Ht  Incident  wave  height,  m  (ft) 

H0  Deepwater  wave  height,  m  (ft) 

Hs  Significant  wave  height 

Hm0  Energy-based,  or  zero-moment,  estimate  of  significant  wave  height,  m  (ft) 


i>M{4>}  Imaginaiy  part  of  4> 

Kd  Stability  coefficient 

Kr  Reflection  coefficient  of  a  solid  boundaiy,  refraction  coefficient 

Kr  coast  Reflection  coefficient  of  a  solid  boundary 

Ks  Shoaling  coefficient 

l  Length  scale 

L  Wavelength 

ljlp  Linear  scale  of  the  model 

La  Deepwater  wave  length 

Lp  Wavelength  for  waves  at  peak  frequency 

m  Model  quantity  (subscript) 

n  Unit  normal  vector  directed  into  the  solid  region 

N  Number  of  units  in  a  given  area 

Na  Normalized  transmission  coefficient 

Nt  Number  of  HARBD  computational  wave  periods  for  spectral  approximation 
p  Prototype  quantity  (subscript) 

r  Radial  polar  coordinate,  m  (ft) 

Re{(j)}  Real  part  of  (j) 
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Cl 


5  Specific  gravity  of  an  individual  armor  unit  relative  to  the  water  in  which  it 
is  placed,  Sa  =  yjyw 

S(fi)  Spectral  energy  density  at  frequency,/ 

Tp  Peak  spectral  wave  period,  sec 

u,v  Horizontal  velocity  components,  m/sec  (ft/sec) 

V  Armor  unit  volume 

Wa  Weight  of  an  individual  armor  unit 

wk  Weighting  factor  for  /fc’th  HARBD  computational  frequency 
x,y  Horizontal  coordinates,  m  (ft) 

P  Dimensionless  bottom  friction  coefficient 

y  Spectral  peak  enhancement  factor;  phase  shift  between  stress  and  flow 
velocity 

ya  Specific  weight  of  an  individual  armor  unit 

yw  Specific  weight  of  water 

Ax  Element  size 

e  Significant  wave  steepness 

6  Wave  phase 

6  Angle  of  structure  slope  measured  from  horizontal  in  degrees 

K  Wave  number 

A  Complex  bottom  friction  factor 

<t>  Velocity  potential;  packing  density  coefficient 

<ff  Velocity  potential  of  the  scattered  wave 

co  Angular  wave  frequency,  radians 

V  Horizontal  gradient  operator 

d  Partial  differentiation  symbol 
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